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NON-LINEAR EIGENVALUE PROBLEMS AND APPLICATIONS 
TO PHOTONIC CRYSTALS 

CHRISTIAN ENGSTROM, HEINZ LANGER, AND CHRISTIANE TRETTER 


Abstract. We establish new analytic and numerical results on a general class 
of rational operator Nevanlinna functions that arise e.g. in modelling photonic 
crystals. The capability of these dielectric nano-structured materials to control 
the flow of light depends on speciflc features of their eigenvalues. Our results 
provide a complete spectral analysis including variational principles and two- 
sided estimates for all eigenvalues along with numerical implementations. They 
even apply to multi-pole Lorentz models of permittivity functions and to the 
eigenvalues between the poles where classical min-max variational principles 
fail completely. In particular, we show that our abstract two-sided eigenvalue 
estimates are optimal and we derive explicit bounds on the band gap above 
a Lorentz pole. A high order flnite element method is used to compute the 
two-sided estimates of a selection of eigenvalues for several concrete Lorentz 
models, e.g. polaritonic materials and multi-pole models. 


1. Introduction 


Many physical systems are passive in the sense that they do not produce en¬ 
ergy. For example, in linear electromagnetic field theory energy can be transferred 
from the electromagnetic field into the material, but not from the material into 
the electromagnetic field. Moreover, materials are in general dispersive which, 
when frequency is the spectral parameter, results in a non-linear operator function. 
Therefore a large number of systems are accurately described by Nevanlinna func¬ 
tions (also called Herglotz functions) whose values are differential operators. In 
numerical analysis matrix-valued Nevanlinna functions such as Schur complements 
are known as versatile tools [PW02] . In recent years operator-valued Nevanlinna 
functions and their analytical properties have been studied intensively [TreOSi Intro¬ 
duction]. However, we still lack a more detailed understanding of a class of rational 
Nevanlinna operator functions that are sufficiently general to cover several im¬ 
portant physical applications, such as acoustic and electromagnetic problems with 
frequency dependent materials e.g. in photonic crystals. Mathematical research in 
this direction has started, but is still in its infancy |EnglO| , [EKE12| . |Schl3) . 

Here we establish, under new and very general conditions, a complete picture 
of the spectral properties for such rational operator functions. The novelty of our 
approach is that it applies, in the case of several poles, to the eigenvalues between 
the poles where classical min-max variational principles fail completely. As a result, 
in applications to photonics, we cover piecewise constant multi-pole Lorentz models 
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periodic on some bounded domain Q = iliU... U JIm, as well as permittivity func¬ 
tions e(-, w) where the linear part of the corresponding operator function has eigen¬ 
values A := below the Lorentz poles oj^ ^ The abstract operator functions we 
consider have the Nevanlinna property in the sense that they are analytic on the 
complex plane, have self-adjoint values on the real axis, and a finite number of poles 
which are real and of first order with non-positive residues, and their derivatives 
are non-negative between the poles. This property enables us to introduce gener¬ 
alized Rayleigh functionals, establish variational principles, and derive two-sided 
estimates for all eigenvalues of this important class of rational operator functions. 

We demonstrate the efficacy of the new theory for unbounded operator functions 
modelling photonic crystals. These dielectric nano-structured materials which are 
used to control and manipulate the flow of light |JJWM08] are commonly mod¬ 
elled by periodic Lorentz permittivity functions (ED with several rational terms. 
Explicit computations show that the abstract two-sided eigenvalue estimates are 
optimal and we derive explicit bounds on the band gap above a Lorentz pole. The 
operator function is discretised with a high order finite element method and several 
concrete examples e.g. for polaritonic materials illustrate the general theory. In 
particular, we compute the two-sided estimates of a selection of eigenvalues and 
we illustrate the accumulation of eigenvalues at the poles and the corresponding 
singular sequence. In most examples a continuous finite element method is used 
to compute the eigenvalues, but in cases were a block diagonal mass matrix is an 
advantage a discontinuous Galerkin method is employed. 

The paper is organized as follows. In Section [5] we set up the required operator 
theoretic framework. In Section [3] we consider the one pole case and establish min- 
max variational characterizations and two-sided estimates for all eigenvalues. In 
Section 0] we generalize the min-max principles to the multi-pole case and identify 
cases where a band gap occurs. In Section Owe apply our abstract results to pho¬ 
tonic crystals with multi-pole Lorentz models (11.11) . Section O contains the nume¬ 
rical finite element analysis for several material models, illustrating different fea¬ 
tures of the abstract results such as the occurrence of an index shift or band gaps. 

Throughout this paper we use the following notations and conventions. All Hil¬ 
bert spaces are separable. For a closed linear operator T in a Hilbert space TL 
we denote by kerT, ranT, p{T), a{T), and crp(T) its kernel, range, resolvent set, 
spectrum, and point spectrum, respectively; the essential spectrum of T is defined 
as (Tess(T) :={A S C : T—X is not Fredholm}. If T is self-adjoint, then AGtTess(T) iff 
A G cr(T) and A is not an isolated eigenvalue of finite multiplicity. Further, for a Borel 
set /CR, we denote by Ci{T) the spectral subspace of T corresponding to the set / 
and, if T is bounded from below and /x<mincress(T), by N{T, p) := dim£(_oo_^] (T) 
the number of eigenvalues of T that are < p counted with multiplicities. 

2. Operator theoretic framework 

Recent applications in nanophotonics require to study the spectral properties 
of operator functions that depend rationally on the spectral parameter. From the 
analytical point of view, there are two possible approaches and often a combination 
of both is most advantageous. 

The first one is to analyze the operator function directly using properties of its 
operator coefficients. Here we consider analytic operator functions S whose values 
are linear operators in a Hilbert space % and which are given by 

(2.1) 5(A) = A - A - B{C - X)-^B*, dom5(A) = dom A CH, A S C \ a{C), 
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where A and C are linear operators in Hilbert spaces "H and "H, respectively, B acts 
from H to H, and B, C are bounded. 

The second approach is to consider a linearization of the operator function S, i.e. 
a linear operator Ain a larger Hilbert space that reflects all the spectral properties 
of S. A particular linearization of S in (|2.ip is the block operator matrix A in the 
product Hilbert space H = given by 


( 2 . 2 ) 


A=(^ 

\B* C) 


dom A = dom A © %. 


In fact, S is the first Schur complement of A] the relations between spectral proper¬ 
ties of S and A summarized in the next proposition are well-known and not difficult 
to check (see e.g. [TreOSl Section 2.3]). Recall that the spectrum of S is defined as 


a{S) := {A e C \ a{C) : 0 e a(5(A))}, 


and analogously for the point spectrum crp(iS) and essential spectrum UessiS) of S. 


Proposition 2.1. Let LL be a Hilbert spaee, H. = with Hilbert spaces H, TL. 

Let A be a elosed linear operator in H and let B : H. ^ H, C : H ^ H. be bounded 
linear operators. Then, for the bloek operator matrix A in (12.21) and the operator 
function S in 

i) cr(A) \ a{C) = cr{S), <7p{A) \ (7(C) = crp(S), and (7ess{A) \ a(C) = (7ess(S); 

ii) if Xq € <7p(A) \ u(C) with eigenvector (uq uo)*, then uq is an eigenveetor 
of S at Aq/ 

iii) if Xq € ap(S) with eigenvector uq, then 
veetor of A at Xq. 

Proof. All claims follow from the Schur-Frobenius factorization (see e.e. ITreOSl 

( 2 . 2 . 12 )]), 

I B{C-X)-^\ fS{X) 0 \ f I 0 

0 / yV 0 C-aJ (^(C-A)-iH* / 

since the outer two factors are bounded and boundedly invertible and so is C — A if 
A y cr(C') (comp. [TreOSl Theorem 2.3.3 ii)]); note that closures in [TreOSl (2.2.12)] 
may be omitted here since B and C are bounded. □ 


A G C\ct(C'), 



Uq 

-(C-Xq)-^B*uq 


IS an eigen- 


The claims of the following proposition were proved in various degrees of gener¬ 
ality (see [ALMS94] . [AL95j . [MS96j . and also [TreOSj l: for the convenience of the 
reader we give a simple proof. 

Proposition 2.2. Let TL be a Hilbert spaee, TL = TL®TL with Hilbert spaces TL, TL. 
Let A be a elosed linear operator in TL with eompaet resolvent, and let B : TL ^ TL, 
C : TL ^ TL be bounded linear operators. Then 

i) the essential speetrum of the block operator matrix A in (1221) is given by 

(2.3) tTess(Al) = CTess(C'). 

If A and C are self-adjoint, then so is A; if, in addition, A is bounded from below, 
then so is A and 

ii) if maxcr)^) < mincr(A), then (max (t(C'), mincr(A)) C p(A); 

iii) if dim TL = <70, then A has a sequence of eigenvalues of finite multiplicities 
aeeumulating only at +oo; 

iv) if c> sup cress(C') and (c, c+(5) C p(C) for some <5 > 0, then (c, c+e) C piA) 
for some e > 0. 
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Proof. The block operator matrix ^ is a bounded symmetric perturbation of the 
block diagonal operator matrix diag(A, C) in Hence A is closed since so 

are A and C, A is self-adjoint if so are A and C, and A is semi-bounded if so is A. 

i) To prove (12.3L we consider the second Schur complement 5 of ^ given by 

S{X) = C - X- B*{A-domS{X)=n, XGC\a{A). 


Since the two outer factors in the corresponding Frobenius-Schur factorization (see 
e.g. [TFiOSl (2.2.11)]) 


^-A 


f I 0\/A-X 0 {A-X)-^B\ 

[b*{A-X)-^ i)[ 0 5(A) Ao I J 


X G C\a(H), 


are bounded and boundedly invertible, it follows that cress(Al) \ cr(A) = cress(5) 
(see |Tre081 Theorem 2.4.7]). Since B is bounded and A has compact resolvent, 
B*(A — X)~^B is compact in H and hence tTess(5) = CTess(C') \ (t{A). Altogether we 
obtain that, for A ^ cr(A), 

(2.4) A G (Tess(-A) > A G fTess(5) 4-> A G 0'ess(C’)- 


It remains to be proved that A G dess (A) 4=4> A G cress(C') also for A G cr{A). In 
this case, since A has compact resolvent, A is an isolated eigenvalue of A and the 
algebraic eigenspace C\{A) of A at A is finite dimensional. If P\ is the orthogonal 
projection in % onto Cx{A) and we choose ^ G M\cr(A), then Aq := A— {X — ^)P\ 
is a finite dimensional perturbation of A with A ^ ct(Ao). Then the equivalence 
(I2.4[l applies to the block operator matrix 


Ao 




dom A = dom A © H, 


and yields that A G cress(Ao) 4=^ A G cress(C')- Since Ao is a finite dimensional 
perturbation of A, we have (Tess(Ao) = dess (A), which completes the proof of (12.31) . 

Now suppose that A and C are self-adjoint, and hence so is A. 

ii) The claim is an immediate consequence of [AL95[ Thm. 2.1] (see also |Tre091 
Theorem 5.2, Corollary 5.3]). 

iii) Since A is semi-bounded with compact resolvent in % and dim H = oo, the 
spectrum of the block diagonal operator matrix diag (A, C) in 77 © "H consists of the 
sequence {iyj{A))JXi of eigenvalues of A of finite multiplicities accumulating only 
at +00 and of the spectrum of C. As A is a bounded perturbation of diag (A, C), 
it will also have a sequence of eigenvalues of finite multiplicities accumulating only 
at +00 (see e.g. [Kat951 Section V.4.3]). 

iv) By the assumptions on A, (7, and c, we know that A has finitely many 
eigenvalues < c and C has finitely many eigenvalues > c (counted with their finite 
multiplicities), say, 


©(A) < i^2(A) < • • • < vn{A) < c < vi{C) < V2{C) < ■ ■ ■ < vm{C). 

If Pj and Qk are the orthogonal projections in 77 onto ker(A — iXj{A)) and in H 
onto ker(C' — i^fc(C')), respectively, and VN+iiA) is the smallest eigenvalue of A that 
is > c, then 


N M 

Ai := A - ^ (i.,(A) - ,xn+i{A))P,, Ci := C - ^ (i/fe(C) - c)Qk 
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are finite dimensional perturbations of A and C with the property that 


maxcr(Ci) = c < iy]v+i(A) = mincr(j4i). 
Now claim ii) applied to the block operator matrix 

Ai := C )’ dom^ = dom A © "H, 


1 


yields that {c,vm+i{A)) C p{Ai). Since the self-adjoint operator A is a finite 
dimensional perturbation of the self-adjoint operator Ai, it follows that A has only a 
finite number of eigenvalues in the interval (c, ujv+i(A)) counted with multiplicities 
(see |BS871 Theorem 9.3.3, p. 215]). Hence (c,c + e) C p[A) for some e > 0. □ 

Remark 2.3. In fact, the total multiplicity of the eigenvalues of A in the interval 
(c, PAr+i(A)) is at most N + M, i.e. at most the sum of the multiplicities of all 
eigenvalues of A less than or equal to c and of those of C greater than c; this 
follows from [BS871 Theorem 9.3.3, p. 215]. 

Since the diagonal part diag (A, C) of A in H H is bounded from below and 
the off-diagonal part is bounded with norm IJiljl, the block operator matrix A is 
bounded from below with 

min{min(T(A), min cr(C')} — jjHjj < min(T(A) < min{min cr(A), min cr(C')} + jjHjj 

(see |Kat95l Theorem V.4.11] and e.g. [CT151 Lemma 5.2]). The off-diagonal nature 
of the perturbation allows one to strengthen this result (see |KMM07] for the case of 
bounded A); the following proposition is immediate from [Tre091 Theorem 5.6 (i)]. 

Proposition 2.4. The block operator matrix A in (12.2|) is bounded from below with 


min{min(T(A), min (t(C')} — S < mincr(A) < min|min(T(A),mincr(C)} 


with 



here 6 < UBj], and 6 < jjHj] if and only if min( t(C') ^ mincr(A). 

In the next two sections, we will study the spectral properties of the block 
operator matrix A in (12.21) and its Schur complement S given by m in greater 
detail. Here we distinguish the cases that S has one pole and more than one pole. 

3. One pole case: variational principles and eigenvalue estimates 

In this section we consider the case that C = i.e. the operator function S 
in the Hilbert space T-L in (|2.1I) has a single pole at the point c G M, 


(3.1) 5(A) = A — A — , dom5(A) = domA c It, AgC\{c}, 

c — A 


where A is a self-adjoint operator with compact resolvent and bounded from below 
and where H is a non-zero bounded linear operator. 
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3.1. Spectral properties of S. To investigate the spectrum of S, we use its close 
relation to the spectrum of the block operator matrix A in (I2.2|l given by 

(3.2) (^omA = domAo^'H 

Since A is self-adjoint and bounded from below, the spectra cr(.4) and (t{S) = 
<7{A) \ {c} are real and bounded from below. By Propositions 12.II and 12.21 we have 


fHess('^) — ^ess(cfi^) — 



U'ess(tS) — (Jess(vd) \ — 


if dim It = oo, 
if dim'H < oo, 

0 , 


and (c, c -I- e) n (t{A) = (c, c -I- e) fl cr(5) = 0 for some e > 0. 

Hence the spectrum of A and of S in the intervals (—oo, c) and (c, -|-oo) is discrete 
and accumulates at most at the right end-points c and -l-oo, respectively. We denote 
the corresponding sequences of eigenvalues, ordered non-decreasingly and counted 
with multiplicities, by C (—oo,c) and {X 2 ,j)^Li C (c,-too), respectively, 

with ni, 712 G No U {oo}: 


u(^) \ {c} = ap(^) \ (cj = a(5) = ap(5) = U (A 2 .,);= i. 

Here ni = oo means that the sequence C [mincr(^),c) is infinite and 

accumulates at c, while 772 = oo means that the sequence (A 2 ,j)Jd]^ C [c-|-£, -|-oo) is 
infinite and accumulates at -|-oo. Since A has compact resolvent, we have 

(3.3) 771 < oo if dim'H < oo, 772 < oo if and only if dim'H < 00 . 

Indeed, if diuv'H < 00 , then cress(Al) = cress(C) = 0 by Proposition 12.21 the second 
claim follows because A, and hence A, is bounded if and only if dim'H < 00 . 


^ 1,1 

-•— 


Ai,2 

— • — 


Ai,3 

— • — 


A2,2 

- • - 


^2,3 

-•- 


Figure 1 . Eigenvalues (Ai,j)N C (—00, c) and (X2,j)Jli C (c,+ 00 ) of S. 


The condition dimH = 00 is only necessary for eigenvalue accumulation at c, 
i.e. for 77i = 00 . The following two propositions show that dim'H = 00 and B 
injective are sufficient conditions for 771 = 00 . Later, with more advanced tools, we 
will be able to show that the weaker condition dim ran H = 00 is both necessary 
and sufficient (see Theorem 13.91 hill. 

Proposition 3.1. Let % he a Hilbert spaee, H = with Hilbert spaces TL, TL. 

Let A be a self-adjoint operator in H with compact resolvent and bounded from 
below, B :% a non-zero bounded linear operator, and c S R. 

If dim H. = 00 and c is not an eigenvalue if infinite multiplicity of A, then c is 
an accumulation point of eigenvalues of A, and hence of S, from the left but not 
from the right. 

Proof. If dimH = 00 , Proposition 12.21 shows that aess{A) = (cj. Since A is self- 
adjoint, this implies that c is either an eigenvalue of infinite multiplicity or an 
accumulation point of eigenvalues. The former is excluded by assumption and 
(c, c -I- e) C p{A) for some £ > 0 by Proposition l2.2l ivl. Thus c is an accumulation 
point of eigenvalues of A, and hence of S by Proposition l2.II iil. from the left but 
not from the right. □ 
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Proposition 3.2. Let % he a Hilbert space, H — T-L^T-L with Hilbert spaces H, H. 
Let A be a self-adjoint operator in H. with compact resolvent and bounded from 
below, B :% ^H. a non-zero bounded linear operator, and c S R. 

If B is injective, then dimker(yt — c) < N(A,c); in particular, c is not an 
eigenvalue of infinite multiplicity of A. 


Proof. First we show that 

(3.4) M= (m m)* G ker(.4 — c) \ {0} => c = ^ 

hr 

In fact, — c)m = 0 is equivalent to 

(3.5) {A — c)u Bu = 0, 

(3.6) B*u = 0. 


If M = 0, then (ixa implies Bu = 0 and thus u = 0 since B is injective by 
assumption, a contradiction to m ^ 0; hence m ^ 0. Taking the scalar product with 
u in (13.51) and using dSH), we find that 


0 = ((A — c)u, u) + {Bu, u) = {{A — c)u, u) + {u, B*u) = ((A — c)u, u), 

which completes the proof of (13.41) . Next we prove that, if P : "H © "H —>■ "H denotes 
the projection onto the first component, then 


(3.7) ric := dimspan(Pker(.A — c)) < N{A,c). 

Assume to the contrary that Uc > N{A,c) + I. We can characterize the eigenval¬ 
ues ui(A) < V 2 {A) < ••• < yN{A,c){^) < c < UAr(A,c)-i-i(^) ^ of A counted 
with multiplicities by the classical min-max variational principle (see e.g. [RS781 
Theorem XIII.I]). If we use Pker(A—c) C dom A and (13.4L we obtain the estimate 


Vj{A) = min max 


{Au, u) 


£Cdoni A 
dim C.= 


¥^0 


,2 ^ 


mm max ■ 

£CP ker(^-c) 

dim C=j u^O 


{Au, u) 


= c, j = l,2,...,nc. 


Since N{A,c) + I < ric by assumption, it follows that V]\{(^a,c)+i{I^) < c, a contra¬ 
diction. 

Since A has compact resolvent, (1X71) implies that Uc < N{A,c) < oo. Thus for 
every choice of ric + I elements Ui = {ui UiY G ker(A — c) \ {0}, i = I, 2,..., ric + 1, 
the first components {ui,U 2 , ■ ■ ■, Un^+i} are linearly dependent. Hence there exists 
{ai)^lh ^ C"'=+^ \ {0} such that aim = 0. Since G ker(A — c), 

(13.41) implies that aim — 0 and hence aim = 0. 

Altogether, it follows that dimker(A — c) < < N{A, c) < oo. □ 


Remark 3.3. i) A necessary condition for B injective is that dimH. < dimP. 
ii) If A > c, then N{A,c) = 0; in this case B injective implies c ^ crp(A). 

Proposition 3.4. If, under the assumptions of Proposition^^ B is injective and 

(3.8) u G dom A n kerP*, u 0 (A —c)u^ranP, 

then c crp(A). Conversely, if c ^ crp(A), then p.8p holds. 

Proof. If cG crp(A), then (13.41) and (13.51) show that the first component of every 
eigenvector at c does not satisfy (j3.811 . Vice versa, if (13.8p does not hold, there 
exists u G dom(A) D kerP*, u 0, and u G H. such that (A — c)u = —Bu. Then 
the non-zero vector {u ft)* is an eigenvector of A at c and hence c G crp(A). □ 
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Remark 3.5. i) If kerB* = {0}, then condition (13.81) is trivially satisfied, 
ii) A sufficient (but not necessary) condition for (13.8p is 

(3.9) M G dom A n ker B*, u 7 ^ 0 {{A — c)u,u) 0; 

in fact, the latter implies that {A — c)u ^ (kerB*)-*- = ranB. 

In the next subsection, we characterize the eigenvalue sequences (Aij))Aj^ and 
(A 2 j )J1i of S by min-max variational principles. The following proposition provides 
necessary information for this and, in particular, for the index shifts occurring 
therein; here it is crucial that the spectrum of S is bounded from below and has a 
gap to the right of c. 

Proposition 3.6. Let H be a Hilbert space, H. — 'HtB'H with Hilbert spaees H, H. 
Let A be a self-adjoint operator in H. with compact resolvent and bounded from 
below, B ^ % a non-zero bounded linear operator, c G R, and 

N{A,c) = dim/:(_oo,c](A). 

Then the operator function S given by (EH) has the following properties: 

i) S satisfies Assumptions (i)-(iv) of |EL04) on (—oo,c) and on (c,+oo); 

ii) there exist oq G (—oo, c), oi G (c, +oo) with (—oo, oq) C p{S), (c, oi) C p{S); 
hi) the dimensions 

dim£(_oo,o)('5(A)) =: Ki, AG(-oo,ao), 

(3.10) dim£(_oo,o)('5(A)) =: K 2 , AG(c,ai), 
are independent of A and finite with 

Ki = 0, K2< N{A, c); 

iv) if B has closed range, P denotes the orthogonal projection onto £(_oo,c](4l) 
and Q the orthogonal projection onto kerB*, then 

(3.11) K 2 < rank (PQ) < min{ N{A,c), dim ker B* }. 

Proof, i) Assumptions (i) and (ii) in |EL04) are satisfied since dom5(A) = dom A is 
independent of A and S is holomorphic on (—oo, c) and on (c, +oo). Because 

(3A2) ^(5(A)u,w) =-lluf - <-||uf, wGdomA, UT^O, 

we have 5'(A) < —I for A G (—oo,c) and for A G (c, +oo). Thus also Assumption 
(iii) in [EL04) is satisfied for S; note that the values —oo and +oo, respectively, 
of the functional p in |EL04[ (2.3)] can be replaced by the end-points a and /3, 
respectively of the interval in which the operator function is considered, here c 
in both cases. For fixed A G K \ {c}, the operator 5(A) is a bounded symmetric 
perturbation of the self-adjoint operator A — A which is bounded from below with 
compact resolvent. This shows that 5(A) has finitely many negative eigenvalues 
and hence Assumption (iv) in |EL04) is satisfied. 

ii) Due to i) we can apply [EL04[ Lemma 2.6] which yields that the function 
A I—>■ dim£(_oo Q) (5(A)) has constant values ki on (—oo,ao) and K 2 on (c, Oi) for 
all ao G (—oo,c), a\ G (c,-l-oo) with (—oo,ao) C p(5), (c, oi) C p(5). 

Since {S{X)u,u) —>■ -l-oo, u G domA, u 0, for A —>■ —oo and 5'(A) < —I for 
A G (— 00 , c), there exists ao G (—oo, c) with 5(A) » 0 for A G (—oo, ao)- Hence we 
have (— oo, ao) C p(5) and ki = dim £(_oo,o) ('5(A)) = 0 for all A G (— oo,ao). 
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Since {c} = cr(C'), Proposition 12.21 ivi shows that (c, c + e) C p{A) \ {c} = p{S) 
for some e > 0, without loss of generality c+e < i^n+i{A), and hence we can choose 
ai = c + e. With P defined as in claim iv) and := I — P, we have 

(3.13) PAP < cP, P^AP-^ > VN+i{^)P^- 
Then, for A G (c, oi) and u € dom^, 

(3.14) (5(A)m,u) = {P{A-X)Pu, u) + {P^{A-X)P^u, u) - 

c — X 

> {P{A-X)Pu,u) 


since c < A < oi < vm+i{A) which implies that the last two terms in (13.1411 
are non-negative. By the classical min-max variational principle, this shows that 
the number of negative eigenvalues of 5(A) is less than or equal to the number of 
negative eigenvalues of P{A — X)P. The latter is N{A, c) by definition of P and so 
K 2 = dim(5(A)) < N[A,c) for A G (c,ai). 

iii) In order to show the first inequality in (I3.11L let := I — Q. We proceed 
from the equality in (13.141) above, but now we keep the last term in (13.141) to obtain 
that, for A G (c, ai) and u G dom A, 

(3.15) {S{X)u,u) > {P{A — X)Pu,u) + -—r-’■— 

A — c 


Since A > c, the operator P(A — X)P is non-positive. Hence, ioi u £ H, u = 
{ui U 2 Y G kerH* 0 ranH, we can estimate 


(3.16) 


> 2 ^[QP{A - X)PQui, ui) + [Q^P{A - X)PQ^U2, ^ 2 )^ 

Since A>ui{A) and A<ai, we have P{A— X)P>{vi(A) — ai)P . Because B has 
closed range by assumption and is the orthogonal projection onto (kerH*)^, 
there exists /3 > 0 such that Q^BB*Q^ > fP on ran 5. Hence, if we choose 
A G (c, oi) such that 

/ 3 " 


(3.17) 


c < A < c + 


then (I3.15L (13.1611 yield that 
> 2 {QP{A - 

>2{QP{A-X)PQui,ui). 


2(ai - vi{A)) ’ 


/32 


{S{X)u,u) > 2(QP(H-A)PQri,ui) +2(pi(H)-ai)||u2f+ ^lk 2 f 

A — c 


Now the classical min-max variational principle shows that, for A with (13.171) . the 
number of negative eigenvalues of 5(A), is less than or equal to the number of 
negative eigenvalues of the negative operator QP{A — X)PQ, which is equal to 
rank(P(3). The second inquality in (13.111) is obvious. □ 


3.2. Variational principles. In the following we characterize the two eigenvalue 
sequences C (—oo,c), (A 2 ,j)jTi C (c,+ 00 ) of 5 by means of min-max 

variational principles. They are based on variational principles for eigenvalues of 
analytic operator functions established in [EL04) . applied to the operator function 
5 in the intervals (—oo,c) and (c,+ 00 ). 

One of the key properties of the Schur complement ensuring that the assumptions 
of |EL04] are met is its Nevanlinna property ensuring monotonicity on (—oo,c) 
and (c,+ 00 ). It guarantees the existence of (at most) one zero of the function 
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A i-A (5(A)u, u) on each of these intervals; this zero serves as a generalized Rayleigh 
functional in the variational principle. 


Lemma 3.7. Let % he a Hilbert space, % = % (B H with Hilbert spaces H, H. 
Let A be a self-adjoint operator in H with compact resolvent and bounded from 
below, B : H ^ H a non-zero bounded linear operator, and c S K. Then, for every 
u G dom^l, u ^ 0, the function 

(3.18) ipu{X) ■■= {S{X)u,u) = {{A - \)u,u) - - - AgK\{c}, 

has at most one pole at c, is strictly decreasing on (—oo, c) and on (c, +oo). Thus (fu 
has at most one zero pi{u) G (—oo, c) and at most one zero P 2 {u) G (c, +oo), given by 


(3.19) 


Pl,2{u) 


1 f {Au, u) 



T 



{Au,u) _ \ 

hP 7 


I7IP 


Remark 3.8. i) If B*u = 0, then is linear with zero {Au, u); nevertheless, formula 
(I3.19|) is still meaningful and gives pi{u) = Timi{{Au, u),c}, p 2 {u) = iiiax.{{Au, u),c}. 

ii) The functionals pi ,2 induced by the numerical range of the Schur complement 
S are related to the functionals A± induced by the quadratic numerical range of 
the block operator matrix A (see |LLT02[ (3.1) and Lemma 3.7]). 

Proof. The claimed monotonicity is immediate from inequality (|3.12|1 . If B*u ^ 0, 
then the claim follows from the fact that the function in (I3.18p has a pole at c 
with {S{X)u,u) —>■ +00 for A —>■ —oo and A \ c, {S{X)u,u) —>■ —oo for A ^ c 
and A —>■ +oo. If B*u = 0, then the assertion is immediate from the fact that the 
function in (13.1811 is linear with precisely one zero {Au,u). 

Formula (|3.19l) follows from the fact that A G K \ {c} is a zero of if and only 
if it is a solution of the quadratic equation 

[{A — X)u, u) (c — A) — {BB*u, u) = 0. □ 

Theorem 3.9. Let H be a Hilbert space, TL = TL (BTL with Hilbert spaces TL, TL. 
Let A be a self-adjoint operator in TL with compact resolvent and bounded from be¬ 
low, B : TL ^ TL a non-zero bounded linear operator, and c G K. Let N{A,c) = 
dim£(_oo c] < dim'H and let the eigenvalues {r'j{A))jTfl^ of A, counted with mul¬ 
tiplicities, be ordered non-deereasingly, i.e. 

{A) < • • • < l'N(A,c) (^) < C < I^Ar(A,c) + l (^) < ' ' ’ • 

Then the spectrum of S consists of two {finite or infinite) eigenvalue sequences 
(Aij)Jli C (—oo,c) and {X 2 ,j){^Ti C (c,+oo), ni, n 2 G No U {oo}, which can be 
characterized as 

(3.20) Ai , = min max pi(m), j = 1,2,..., m, 

’ £Cdom A u^C 

dim £,=j u^Q 

(3.21) ^ 2 , 1 = min max p 2 (u), j = 1, 2,..., n 2 . 

Z^CdomA u^C- 
dim C=j + K2 u^O 

Here pi, p 2 are the functionals in (13.191) and the index shift K 2 defined in (13.101) 
satisfies K 2 < N{A,c); in partieular, K 2 = 0 if A > c. Moreover, 

i) if B has closed range and P and Q are the orthogonal projections onto 
B(-oo,c]{A) and^^e^:B* respectively, then 

K 2 < rank (PQ) < min{A^(24, c),dimkeri3*}; 
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ii) 712 = 00 if and only if dimH = c»; in this case, (A 2 j)^i C (c,-l-cxj) accumu¬ 
lates {only) to + 00 ; 

iii) ni = 00 if and only if dimrani? = oo; in this case, {^ij))^i C {—co,c) 
accumulates {only) to c from the left. 

Proof. Proposition 12.11 yields that crp(5) = crp{A) \ {c}. Moreover, ni + 712 + 
dim£{c}(-4) = dim'H + dim'H. 

By Proposition 13.61 ii. the Schur complement S satisfies Assumptions (i) to (iv) 
of [EL041 Theorem 2.1] on each of the two intervals (—oo,c) and (c, + 00 ). Hence 
the variational characterizations (|3.20p and (I3.21|) follow from |EL04[ (2.9] applied 
on (—oo,c) and (c,+ 00 ), respectively; the claims for the index shifts ni in (—oo,c) 
and K 2 in (c,+ 00 ), in particular i), follow from Proposition 13.61 iii. 

Claim ii) was proved in (13.311 . For the proof of claim iii), we first suppose that 
dim ran B < 00 . Then ^ is a finite rank perturbation of diag {A, c) and hence there 
are at most finitely many eigenvalues in (—oo,c), i.e. ni < 00 . For the converse, 
assume that dimrani3 = oo. If a; ^ kerH* = ker{BB*), then 

(3.22) lim (5'(A)a;, x) = — 00 . 

A/'c 

Now let j S N be arbitrary. Then there exists a subspace £j C dom A with dim £j — 
j and £j Piker = {0}. In fact, since ranH is infinite dimensional by assumption, 
there exists a subspace £' C ranH with dim£' = j. Let {e'^, 62 ,..., e'} be an 
orthonormal basis of £' and e G ( 0 , 1 / 2 ). As dom A is dense in H, there exist 
{ei, 62 ,..., Cj} C dom A such that \\ek — e/|| < 6 , k = 1,2,... ,j. Then 

(3.23) |(ei,e/)| < +e, i,k = 1,2,... ,j, 

where Sik is the Kronecker symbol. Since e < 1/2, it is easy to see that (I3.23|) implies 
that {ei, 62 ,..., Bj} are linearly independent and that £j := span {ci, 62 ,..., e^} ("I 
kerH* = spanjei, 62 ,..., e^} n (ranH)-*- = { 0 }. 

Since £j n kerH* = {0} and £j is finite dimensional, whence closed, relation 
(I3.22|) implies that there exists Aq < c such that (5(Ao)x, x) < 0, x G £j \ {0}. The 
latter means that 

Pi(x)<Ao, xG£j\{0 }. 

It follows that max 3 .g£^\{o} Pi{x) < Aq. Now the variational principle (13.201) implies 
that Xj < Aq. Therefore there are at least n eigenvalues of S below c. Since n was 
arbitrary, there must be infinitely many eigenvalues of S below c. □ 

Remark 3.10. The variational principles (13.201) . (13.211) continue to hold if we replace 
the domain dom A of the operator A therein by the form domain dom a = dom j Aj^/^ 
of the quadratic form a associated with A, i.e. a[M] := {Au,u) (see [KLT041 Theo¬ 
rem 3.1, Lemma 3.5]); note that the functionals Pi, 2 (m) in (13.191) are defined for all 
u G dom a. 

Remark 3.11. Propositions 13.11 and 13.21 together imply that if dimH = 00 and B 
is injective, then c is an accumulation point of eigenvalues of A, and hence of S, 
from the left. Theorem 13.91 iiil shows that dim ran H = 00 ensures eigenvalue 
accumulation at c from the left. 

The latter result is stronger because dim'H = 00 and B injective imply that 
dimranil = oo. Otherwise, if dimrani3<oo, there exist ns G N and biGR, biGR, 
i = l,2,... ,nB, with B = YfJi^i{-, bi)bi. Since dim'H = 00 , there exists uo€'H\{0}, 
Uq J- bi, i = l,2,..., ns. Then Buq = 0, a contradiction to B being injective. 
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The estimate for the index shift K 2 in Theorem 13.91 i) is of special interest for 
numerical approximations where infinite dimensional spaces have to be replaced by 
finite dimensional ones (see Section 0). 

Remark 3.12. liTL-TL are finite dimensional, then dimker il* = dim "H—dim rani?* = 
dimii — dim(keri?)-*- > diui'H — diinii. Hence the condition dimii — dimT-L > 
N{A,c) ensures that min{iV(H, c), dimkeri?*} = N{A,c). 

So to obtain correct information on the index shift K 2 of an infinite dimensional 
problem with dimkerH* = oo using finite-dimensional approximations in 

one should ensure that dimi^M —dimiiM > N{Am, c) where N{Am, c) 
is the number of eigenvalues of the corresponding operator Am in 'Hm that are < c. 


3.3. Two-sided eigenvalue estimates. Next we show that the min-max varia¬ 
tional principles in Theorem 13.91 provide two-sided estimates for all the eigenvalues 
of S. The bounds are expressed in terms of the eigenvalues of the left upper entry 
A of A and the norm of the off-diagonal entry B. 


Theorem 3.13. Let % he a Hilbert space, % = with Hilbert spaces H, H. 

Let A he a self-adjoint operator in H with compact resolvent and bounded from be¬ 
low, B : 'H ^ H a non-zero hounded linear operator, and c G K. Let N{A,c) = 
dim£(_oo,c] < dimlt and let the eigenvalues of A, counted with mul¬ 

tiplicities, he ordered non-decreasingly, i.e. 


( 3 . 24 ) i^l(kl) < ■ • ■ < ’^N(A,c){^) < C < I^Ar(A,c) + l(^) < ' ' ’ • 

Then the eigenvalues (Aij)",ij^ C (—oo, c), i^2,j)'j=i'^ic,-\-oo), ni, n2 GNoU{oo}, of 
S satisfy the two-sided estimates 

( 3 - 25 ) Xi j < Xij < XYj, X2J+K.2 — 

for j = 1 , 2 ,.. ., ni and j = 1 , 2 , ..., 112 — K2, respectively, where 


(3.26) At,- := 




(3.27) At, := 


AA) 


and 

(3.28)At 


J + K2’ 2 



f ^3+^2iA)-C 


+ \\BP, 


-I-miner (Hi?*) < min{ie, (H), c} 


2 

-I-mintT(i?i?*) > max{H',-+^^(H), c}, 


(3.29)A; 


2,i + K2’ 


Vj-\-K2iA) + c f i/,+«2(^)-c 




Here K 2 and ni, n 2 satisfy 

i) K 2 < N, and K 2 < min{N^(H, c), dimker i?*} if B has closed range; 

ii) n 2 = 00 if and only if dim H = 00 ; 
hi) ni = 00 if and only if dim rani? = 00 ; 
iv) 711-1-712-1- dim£|c}(.4) = dimii -1- dimii, 

m > N{A, c), n 2 > dim H — N{A, c). 


Proof. In order to prove (13.251) . we use the variational characterizations (I3.20L 
(13.211) of Aij, X 2 J proved in Theorem 13.131 with pi _2 as defined in (I3.19|) : note that 
for (Ai,j)”i^, the index shift is 0 by Theorem 13.131 
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For estimating pi ^2 we note that the functions 


(3.30) 


s - 1 - c 

2 +V 


s G M, jS G [ 0 , H-cxd) 

(3.31) 

g{s,P) 

s - 1 - c 

■“2 V 


s G M, jS G [ 0 , H-cxd) 


are both increasing in s; moreover, / is increasing in /3, while g is decreasing in f3. 

These monotonicity properties, together with the classical min-max characteri¬ 
zation of the eigenvalues Vj {A) and the two-sided bounds 


0 < Tmna{BB*) < < \\B*f = \\Bf, 


u^O, 


yield the claimed estimates (13.251) ; note that the leftmost bound yields the estimates 
in (I3.27p . (j3.28l) since e.g. 


\U ^ (^) + c 

^ - 2 - 



= min{uj(A), c}. 


Claim i) follows from Proposition 13.61 Claims ii), iii) were proved in Theorem 
13.91 hi. iii). The first claim in iv) is obvious; the lower bounds for ni and n 2 follow 
from the upper estimate in (13.3511 and the lower estimate in (13.381) . □ 


The following corollary shows how to obtain two-sided computable bounds for 
the eigenvalues and (A 2 ,j)”di of A. It is an immediate consequence of 

Theorem 13.131 and of the monotonicity of the functions f,gin (13.301) . (13.311) . 


Corollary 3.14. Suppose that, for some N < dirnTf, the first N eigenvalues of A 
admit two-sided estimates 

(3.32) i^fiA)<u,iA)<ufiA), j = 1,2,..., iV. 

Then 

(3.33) Xi j < Xij < XYp < X 2 J < , j = 1,2,..., iV, 

where Af^, X^j are obtained from Af^, X^j by replacing Vj{A) by its lower bound 
ly^iA), while XY^j, X^j are obtained from XYj , XYj by replacing Vj [A) by its upper 
bound (^) ■ 

Remark 3.15. If the two-sided bounds for i^jiA) in p.32l) are computable, then 
so are the two-sided bounds for Ai.^, X 2 J in (13.331) . Computable upper bounds 
for Vj{A) may be obtained by the classical min-max variational principle (see e.g. 
[KS84) 1. e.g. from Galerkin approximations since, for every TV-dimensional subspace 
T-Lff oi dom A C H, 

vAA) = min max (Au.u) < min max (Au.u) =: vY (A), 7 = 1, 2,..., TV. 

dimC=j 


Proposition 3.16. Let the assumptions of rTieorem l3.13l hold, set I'oiA) := — 00 , 
and 


(3.34) 


SbA) ■■= ||B|| tan 


1 2||B|| 

- arctan --; 

2 \t-c\ 


t G (— 00 , -foo). 
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Then the eigenvalues C (—c», c) of A satisfy the estimates 

(3.35) Vj{A) - 6 b{vj{A)) < \ij <Vj{A), j = 1,2,... ,N{A,c), 

(3.36) c-5 b(vj(,A)) <\ij <c, j = N(A,c) + l,...,ni, 

and the eigenvalues {^ 2 ,j)^LiC{vi{A),+oo) of A satisfy 

(3.37) c< X 2 ,j <c + SBivj+K 2 iA)), j = l,2,...,N(A,c) - K 2 , 

(3.38) 

Vj+f,^{A) < \2j < i/j+„2(A) + (5B(^’j+K2(^))> j=N{A,c)-K2 + l,...,n2-K2. 


Proof. The upper bound for Xij follows from the definition of X^j in (13.271) : the 
lower bound for X 2 J follows from the definition of in (13.281) and the estimate 

mina-{BB*) > 0 which yield 


> max{i/j+„,(A),c} 


c, j = l,2,...,N{A,c)- K 2 , 

Vj+f,^{A), j = N{A,c) - K 2 + l,...,n 2 - K 2 - 


The lower bound for Aij- and the upper bound for X 2 J follow from the following 
alternative formulas for the solutions of quadratic equations on the right hand sides 
of (13.261) . (13.291) (see e.g. |Tre091 Lemma 5.1 and (5.1)], 

(3.39) Xf j =mm{i'j{A),c} - 5B{iyj{A)), j = 1,2,... ,ni, 

(3.40) >^2,j+K2 = max{j/j+,,2(yl),c} + ^B(^'i+K2(-4)), j = 1,2,... ,n2 - K 2 , 

together with the definition of N{A,c) (see (13.241) 1. □ 


Remark 3.17. Note that, if ni = 00 and 71-2 = 00 , respectively, (13.361) and p.38l) 
also yield the order of convergence of Xij yAc and of X 2 J — i'j+K 2 (^) \0 for j —>■ 00 , 

c-Ai., “0(i.,(^)-i), j>N{A,c) + l, 

X2,j-1'j+^,2{^) "=°°C )(^'4 + K2(^)“^): j>N{A,c)-K2 + l, 

since a Taylor series expansion shows that (5 b( f) = 0{t~^), t —> 00 . 

Classical perturbation theory for the block operator matrix A in (|3.2I) with di¬ 
agonal entries A, c and off-diagonal entries B, B* shows that 

(3.41) ^ a(Al)nK(24,e)(Al) + ||i?||,c-||i?||)=0, 

(3.42) IIBII < ^ (T(Al)n(c+||i3||,^^(24.c)+i(4l)-||i?||)=0. 

However, the off-diagonal structure of the perturbation allows for stronger inclu¬ 
sions. More precisely, the two-sided estimates in Theorem l3.13l and Proposition l3.16l 
provide tighter estimates for subintervals of (i.'jv(^_c)(4l), c) and (c, J^iv(A,c)+i(^)) to 
be free of eigenvalues. 

Proposition 3.18. Let TL be a Hilbert space, TL = TL (B TL with Hilbert spaces 
H, TL. Let A be a self-adjoint operator in TL with compact resolvent and bounded 
from below, B : TL ^ TL a non-zero bounded linear operator, and c G R. If 

N := N{A,c) = dim£(_oo,c](^) 
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and the eigenvalues of A, counted with multiplicities, are ordered non- 

decreasingly [see (I3.24|l ). then the eigenvalues C (—oo,c) and C 

(c, +oo), ni, 712 S No U {oo}, of S satisfy the inclusions 

(3.43) [Xi,j)f^i C {ni{A)-SB[vi[A)),VN[A)), (c-^b(un-i-i(^)), c), 

(3.44) (A 2 j)^“i''^ C (c, c+dB[i^N[A ))), [X 2 ,j)]lN-K 2 +i {vn+i[A),oo) . 

In particular, i/||i3|| < {c—VN{A))'^ + [c—nj\j{A))[c—VN+i[A^))j then 

(3.45) (7{A) n {vn{A), c - 6b[vn+i[A))) = 0, 

and if ||B|| < ^/{c-i 2 n+i[A))^+ {c -VNiA)){c-UN+i[A)), then 

(3.46) a[A) n (c + Sb[vn[A)), vn+i{A)) = %■, 

7 /||i3|| < •\/2niin{c — UAr(^), ujv_|_i(A) — c}, then both (I3.45P and (13.461) hold. 

Note that the lower bound min{ui(A), c} — 5b[vi{A)) for the whole spectrum of 
A in (I3.43|) coincides with the bound in Proposition 12.41 since therein mincr(Al) = 
viiA), min(T(C') = c, and 6 = Sb[vi[A)). 

The eigenvalue estimates (13.451) , (13.461) are better than the classical perturbation 
results (13.411) . (13.421) for two reasons; firstly, they show that the spectral gaps 
[vj[A),c), [c,vn+i{A)) between the diagonal elements A and c may shrink only 
from one side and, secondly, the norm constraint on the perturbation is not only 
improved from a factor A to 1 as a result of the one-sided shrinking, but even to \/2. 

Proof of Proposition 13.181 The first inclusion in (13.431) follows from (13.351) , the fact 
that the function t ^ t — 5B[t) is increasing on (— oo, -koo), and i^i[A) < Vj[A) < 
J^n[A) < c, j = 1,2,...,TV. The second inclusion in (13.431) follows from (I3.36L 
the fact that the function t i-A- 6B{t) is decreasing on (c, -|-oo), and c < un+i[A) < 

Vj{A), j = N \,N + 2, _The inclusions in (13.441) follow from (I3.37|l . (13.381) by 

similar arguments. 

It is not difficult to check that 6b[i^n+i[A)) < c — UAr(^) if and only if ||i?|| < 
c—iyN{A))'^ + {c—i 2 N{A)){c—iyN+i{A)). Hence the latter condition guarantees 
that the two intervals in ()3.43ll are disjoint and thus (13.451) follows. In a similar way, 
claim p.46p follows from (13.441) . Finally, for the last claim it suffices to note that 
e.g. a/2 min{c-Uiv(^), ujv+i iA)-c} < ^J (c- un(^))^ + (c- UAr(H))(c- uat+i (H)). □ 

Remark 3.19. If c £ crp(H), we have c = vn{A) since N = N{A,c) and hence 
Sb{>^n{A)) = ||H|| limi-_>oo tan (A arctan t) = ||i3||. In this case, (13.431) is not appli¬ 
cable and (13.341) becomes 

a{A) n (c-k \\B\\,i 2 n+i{A)) = 0. 

As we have seen above, many of the previous results can be strengthened if 
A > c. In this case no eigenvalue of A lies below c and hence N = N(A, c) = 0, 
which implies that also K 2 (< N[A,c)) = 0 for the index shift K 2 - The following 
remark summarizes these stronger results. The proposition below adds a so-called 
half range basis result for this particular case. 

Remark 3.20. Suppose that, in addition to the assumptions of Theorem 13 .1 31 fand 
hence of Proposition 13.181) . we have A > c. Then the point spectrum of A consists 
of two sequences of eigenvalues of finite multiplicities, (Ai C (—oo, c), tending 
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to c if ni = oo, and C (i/i(A),+oo), tending to +c» if n2 = oo, which 

satisfy the two-sided estimates (13.251) with iV = K2 = 0 , in particular, 

a-{A) n (c, vi[A)) = 0. 

In the following proposition we formulate basis properties for the first compo¬ 
nents of the eigenvectors of A. Observe that, for eigenvalues different from c, the 
latter coincide with the eigenvectors of S. 

Proposition 3.21. Let % he a Hilbert space, H = H (B LL with Hilbert spaces 

H, H. Let A be a self-adjoint operator in H with compact resolvent and bounded 

from below, the sequence of eigenvalues of A ordered non-decreasingly, 

B : H ^ H a non-zero bounded linear operator, and c G M such that A > c. Then 
each of the following sets can be chosen to form a Riesz basis in H: 

i) the first components of the eigenvectors of A corresponding to eigenvalues 
in (—oo,c], i.e. corresponding to {Xi,jYjLi and to c if c & crp(^); 

ii) the first components of the eigenvectors of A corresponding to eigenvalues 
in (c,-foo), i.e. corresponding to 

Proof. The second claim follows directly from |AL951 Theorem 3.5]. The first claim 
follows from |AL951 Theorem 2.3 (iv)], which shows that the restriction A- of the 
block operator matrix A to the spectral subspace corresponding to the spectrum 
in (—oo, c) is unitarily equivalent to a self-adjoint operator in some Hilbert space, 
and from the fact that cr{A-) = cr{A) ("I (—oo,c) = is discrete. □ 


4. Multi-pole case: variational principles 

In this section we consider the case that the rational operator function S in the 
Hilbert space TL in (12.ip is of the form 

(4.1) S{\) = A—A — BiBg ^ dom5(A) = dom A, A G C\{ci, C 2 ,..., cl}, 

« C£ — A 

with Cl, C 2 ,..., CL G R, Cl < C 2 < • • • < Cl, and non-zero bounded linear operators 
Bi from Hilbert spaces Hi to H, t = 1, 2,..., L. As before, A is supposed to be 
self-adjoint and bounded from below with compact resolvent. 

In the sequel, we set cq := —oo, cl+i := -foo and introduce the bounded operators 

Te:= BeB; l = l,2,...,L. 


4.1. Spectral properties of 5. Let'H:='Hi©'H 2 ©-•-©"Hl. Then 5 is the Schur 
complement of the operator matrix A in (12.21) in'H—'HtB'H given by 


(4.2) 



Bl 

B 2 

•••Sl\ 

Bl 

Cl 

0 

... 0 

B*2 

0 

C2 

... 0 

\Bl 

0 

0 

••• CL J 


dom A = dom A © TL, 


i.e. B=[Bi B2 ... Bl) :'Hi®'H2®- ■ •©'Hl-TH, C' = diag (ci/.^^, 02 /.^^,... ,CLlpiJ. 

Since A is self-adjoint and bounded from below, the spectra cr(A) and cr{S) = 
<7{A) \ {ci, C 2 ,..., Cl} are real and bounded from below. By Propositions 12.11 
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and 12.21 we have 

(4.3) o-ess(.4) = cress(C) = {c£ : dim 'Hg = oo, i = 1,2,... ,L}, 

aess{S) = 0-ess(.4) \ Cr(C') = 0, 

and {cl,Cl + £l) C p{A) for some £l > 0. 

Hence the spectrum of S outside of its poles cg, i.e. in the intervals (— 00 , Ci), 
{cg-i, cg), £ = 2,3,..., L, and (cl, + 00 ), is discrete and accumulates at most at the 
points ci,C 2 , ■ ■ ■ ,cl and at + 00 , respectively. Moreover, since A is bounded from 
below and maxcress(C) = cl, Proposition 12.21 shows that eigenvalues of A cannot 
accumulate from the right at the points cq = —00 and cl- 

Unlike the one pole case, we now also have to study the question of eigenvalue 
accumulation between two poles. Our first result shows that there are gaps in the 
spectrum to the right of all points cg, £ = 0,1,..., L, and hence the eigenvalues in 
each interval {cg-i, cg) cannot accumulate at cg-i from the right, £ = 1,2,..., L. 


Proposition 4.1. Let % be a Hilbert space, A a self-adjoint operator in %, bounded 
from below and with compact resolvent, let Tg, £ = 1,2, ..., L, be non-zero bounded 
self-adjoint non-negative operators inH, ci,C 2 , ■ ■ ■ ,cl € K with ci < C 2 < ■ ■ ■ < cl, 
and set cq = — 00 , cl+i = + 00 . Then the operator function S given by 


(4.4) 5(A) = H —A —dom5(A) = domH, A G C \ {ci, C 2 ,..., cl}, 

. Cg — X 
t=i 


has the following properties: 

i) S satisfies Assumptions (i)-(iv) of |EL04] on (cg-i,cg), £ = 1,2,... ,U + 1; 

ii) there exist ag-i € {cg-i,cg) sueh that (c^_i, a^_i) C p{S), £ = 1,2,..., L+1; 
further, 

dim£(_oo,o)('5(A)) =: Ki, AG(-oo,ao), 

dim£(_oo^o)('5(A)) =: Kf, XG{cg-i,ag-i), i = 2,3,..., L-\-1, 


are independent of X and finite with 

^ r 

Ki = 0, Kg<N(A^^ - i -£ = 2, 3,..., L-I-1, 

in particular, kl-li < N(A,cl); 


hi) ifVg-i has closed range for some £ € {2,3,..., L}, Pg-i denotes the or¬ 
thogonal projection onto £(_oo,c£_i] (^) o.nd Qg-i the orthogonal projection 
onto kerP^-i, then 


(4.5) 


Kg < N 
in particular, 


{Qt-i{pg-iAPg-i — ■ 


Cg-i 


(4.6) kl+1 < rank(PLQL) < min {iV(H, cl), dimkerPL}. 


Proof, i) The proof that S satisfies Assumptions (i)-(iv) of [EL04] on each (c^_i, c^), 
£ = 1,2,...,L + 1, is completely analogous to the proof of Proposition 13.61 ib Here 
we note that, since P^, £ = 1, 2 ,..., L, is non-negative and bounded, we still have 
S'{X) < —I and the operator 5(A) is a bounded symmetric perturbation of A — A 
for A G {cg-i,cg). _ 

ii) Due to i) we can apply |EL041 Theorem 2.1 and Lemma 2.6] which yield 
that there are ag-i G {cg-i,cg) with {cg-i,ag-i) C p(5) and that the function 
A I—!► dim£(_Qo_o) (‘^(A)) has a constant value Kg on (cf_i, (Xg-i), £ = 1,2,...,L-|-1. 











18 


CHRISTIAN ENGSTROM, HEINZ HANGER, AND CHRISTIANE TRETTER 


The proof of ki = 0 is the same as in the one pole case since we again have 
(5 (A)m, u) — >■ +00, u G dom A, m ^ 0, for A— >■ —oo, S'{X) < —I for Ag (—oo, ci), and 
so 5(A) 0 for Ag (—oo, ao)- 

In order to prove the upper bound for ^ = 2,3 ..., L + 1, let A G (c^-i, Ci). 
Then the estimate 

(4.7) 5(A)=Al-A-^^^>Al-A-^^^=:ro(A) 

shows that the number of negative eigenvalues of 5(A) is less than or equal to 
the number of negative eigenvalues of To (A). Choosing A = C£_i + e with e > 0 
arbitrarily small, we find 


(4.8) 


L 

To(c£_i +e) = ^ —C£_i— 

j=i 






If we let e \ 0, the bound for ki claimed in ii) follows. 

hi) Let QgLi ■= I — Qi-i be the orthogonal projection onto ranr^_i. Without 
loss of generality, we may assume that eg-i > 0 is so small that (c;_i, Ci_i +et-i] fl 
<j{A) = 0. Then, for A G (cf_i,c^_i + e^-i), we have the estimate 


(4.9) 


5(A) > P,_i(Al - A)P,_i - ^ 


Ti-1 
A — ci-i 


=:ri(A) + 


r^-i 
A — ct-i 


Since ct-i < X < ce, P£-i{A — A)Pi_i is non-positive while Tj/(cj — A) is non¬ 
negative, j = £,£ + 1,..., L. Thus 7i(A) is non-positive and hence, for u G 'H, 
u = (mi U 2 )* G kerr|_j^ © ranT^-i, we can estimate 

{r.(AK„) = (( <3'-T.(A)0tA ('“iW-O'j 

(4.10) VV(Q^-i'^i(^)Q^-i) QtiTlWQtJ \U 2 JJ 

> 2^((3f_iTl(A)(5£_iUi,Ui) + (Q^_iTi(A)(5^_iM2,'U2))- 
Since A > i^i(A), A < C£-i + e^-i, and = 0 if i'i{A) > c£-i + e^-i, we have 


Pe-i{A - X)P£-i > min{(zAi(A) - (cf_i + e^-i)), 


L 


E 


j=e 



l|r,-|| 

ic£-l + S£-i) 


Therefore, for U2 = Q\_£U2 G ranT^-i, 

{Ti{\)u2,U2) 

(4.II) > -max{c^_i+e^_i-i/i(A4),0}||P^_iM2r-^--rlk2|P 

Cj-{c£-i + e£-i) 

> -(max{c^_i+e^_i-i/i(A),0} + y^-- r)lk2f =: - 7 IIM 2 IP 


where 7 > 0. Because T^-i has closed range by assumption and Qj^i is the ortho¬ 
gonal projection onto (kerT^-i)-*-, there exists 7^-1 > 0 with > 7 ^-i 

on ranr^_i. Altogether, ()4.9I1 . (I4.10L and (I4.11L yield that 

{S{X)u,u) > 2(Q^_iri(A)Q,_iMi,wi) - 27 ||w 2 f + E""' 

A - C£_i 
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Hence, if we choose A G (c^_i, c^_i + ££-i) such that 
(4.12) c^_i < A < c^_i + ^ if 7 > 0, 

we have the estimate 

(S(A)u,u) > 2(Q£_iTi(A)Qe-iUi,ui). 

Now the classical min-max variational principle for semi-bounded operators shows 
that, for A as in (I4.12|) . the number K£ of negative eigenvalues of 5(A) is less than 
or equal to the number of negative eigenvalues of Qf-iTL(A)(5^_i. Since ki is 
independent of A, we may choose A = ci-i + e with e > 0 arbitrarily small and 
proceed in the same way as in (14.81) to see that, in the limit e 0, the number of 
negative eigenvalues of Q£-i is given by the right hand side of (14.51) . 

The last claim for £ = L -|- 1 follows since in this case 71 (A) = A — X and the 
operator QlPl{A — X)PlQl is negative for X > cl and thus has rank(Pi(5i) 
negative eigenvalues. □ 

Remark 4.2. If there is just one pole, i.e. L = 1, then the claims in Proposition 14.II 
coincide with those in Proposition 13.61 note that kerP^ = keri?|. 


Af-1,1 • • • Xt^i 

H • ••H •- 

Cl . . . C£_1 



I 

•Cl 


Figure 2. Eigenvalues (Xej)"!^ C {ce-i, ce) of 5. 


According to Proposition l4.11 the spectrum of A and of 5 in all intervals (c^_i, c^) 
may accumulate at most at the right end-points for 7 = 1,2,..., L-l-1. We denote 
the corresponding sequences of eigenvalues, ordered non-decreasingly and counted 
with multiplicities, by C (c^_i,Cf) with G Nq U {oo}: 

L + l 

a{A) \ {ci,C 2 ,... ,cl} = crp(^) \ {ci,C 2 , ...,cl} = a{S) = ap{S) = |J (X£j)ji^. 

e=i 

Here n£ = oo means that the sequence (Xej)^^i C (c^_i,c^) is infinite and accu¬ 
mulates at C£. It is not difficult to see that 

ne < oo if dim Ri < oo, £ = 1,2,... ,L, 
riL+i < oo if and only if dim'H < oo. 

Both claims are consequences of the fact that A has compact resolvent. Indeed, 
if AivaT-Li < oo, then ci ^ cress(^) by Proposition 12.21 the second claim follows 
because A, and hence A, is bounded if and only if dimH < oo. 

The condition dim 7^^ = oo is only necessary for eigenvalue accumulation at C£, 
i.e. for n£ = oo. The following two propositions show that dim77^ = oo and 
Bi injective are sufficient conditions for = oo, in analogy to Propositions 13.11 
and l3.2l where L = I. In Theorem ld.Gl iiil below, we show that the weaker condition 
dim ran 77^ = dim ran Pf = oo is both necessary and sufficient. 

Proposition 4.3. Let % he a Hilbert space, H = H (B Hi (B 'H 2 (B ■■■ (B Hl with 
Hilbert spaces H, Hi, H 2 , ■■., Hl- Let A be a self-adjoint operator in H with 
compact resolvent and bounded from below, Bi : Hi ^ H, i = 1,2,..., L, bounded 
linear operators, and ci,C 2 ,... ,cl G M with ci < C 2 < ■ ■ ■ < cl. 

If dim77^ = 00 for £ G {I,2,...,L} and c is not an eigenvalue of infinite 
multiplicity of A, then ci is an accumulation point of eigenvalues of A, and hence 
of S, from the left but not from the right. 
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Proof. The proof of Proposition 14.31 is analogous to the proof of Proposition 13.11 if 
one replaces Proposition 12.21 ivl by Proposition 14.II iil. □ 

Proposition 4.4. Let % be a Hilbert space, H — H (B Hi (B 'H 2 (B ■ ■ ■ (B Hl with 
Hilbert spaces H, Hi, H 2 , ■■■, Hl- Let A be a self-adjoint operator in H with 
compact resolvent and bounded from below, B( : He —>■ H, £ = 1,2,..., L, bounded 
linear operators, and Ci,C 2 , ■ ■ ■ ,cl £ K with ci < C 2 < ■ ■ ■ < cl■ 

If Be is injective for £ £ {1,2 ,...,L}, then diniker(^ — q) < N{A,ce); in 
particular, ce is not an eigenvalue of infinite multiplicity of A. 

Proof. The proof of Proposition HU] is completely analogous to the proof of Propo¬ 
sition [321 and thus left to the reader. □ 

Proposition 4.5. If, under the assumptions of Provosition 13131 Be is injective for 
£ € {1, 2 ,..., L} and 

(4.14) u G domAn ker B( , u ^ 0 => {A — ce)u ^ van (^Bi B 2 ... Bl) ■ 
then Ce ^ tTp(.4). Conversely, if ce ^ a-p{A), then ()4.14l) holds. 

Proof. The proof of Proposition l4.5l is completely analogous to the proof of Propo- 
sition [321 and thus left to the reader. □ 

4.2. Variational principles. In this subsection we establish min-max variational 
principles for the eigenvalues of the operator function S in (14.41) in the intervals 
(ce-i,ce), £=l,2...,L + l. 

Theorem 4.6. Let H be a Hilbert space, A a self-adjoint operator in H, bounded 
from below and with compact resolvent, let Fe, £ = 1,2,..., L, be non-zero bounded 
self-adjoint non-negative operators inH, ci,C 2 ,... ,cl £ K with ci < C 2 < ■ ■ ■ < cl, 
and set cq = — 00 , cl+i = -foo. 

Then the operator function S in (HU) is strictly decreasing in ]R\ {ci, C 2 ,..., cl} 
and hence, for each u £ dom(A) = dom(iS(A)), the function X i-A- (S{X)u,u) has at 
most one zero in each interval (c£-i, ce), £ = 1, 2,..., L -|- 1. If we define pe{u) £ 
[ce-i, Ce] for u £ dom(A) = dom(iS(A)) by 

{ Xe{u) if {S{Xeiu))u,u) = 0 for Xe{u) G {ce-i,ce), 
ce-i if {S{X)u,u) <0 for all X G {ce-i,ci), 

ce if {S{X)u,u)>0 for all XG(ce-i,ce), 

then the spectrum of S consists of L -\- 1 {finite or infinite) eigenvalue sequences 
{Xej){l‘Li C {ce-i,ce), £ No U { 00 }, which may be characterized as 

(4.16) Xe j = min max pf(it), j = 1,2,... ,ne, £ = 1,2,..., L-\-1. 

’ Z^Cdom A u^C 

dim C=j + K£ u^O 

Here the index shifts Ke defined as in Proposition 14.11 are finite and satisfy the 
estimates therein; in particular, ki = 0 and kl+i < N{A,cl). Moreover, 

i) if LL has closed range, Pl is the orthogonal projection onto £(_oo,c_l] (^) 
and Ql the orthogonal projection onto kerP/,, then 

kl +1 < rank{PLQL) < min {iV(A, cl), dimkerPi}; 

ii) nL+i=oo if and only if dim'H = oo; in this case, (Ai,+i,j)^i C (c^,-|-oo) 
accumulates {only) to -l-oo; 

iii) ne = 00 if and only if dimranP^ = 00 , £ = 1,2,..., L; in this case, 
(Af,j)^i C (cf_i, c^) accumulates {only) to ce from the left. 
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Proof. According to ProDOsition l2.1l ii. we have cr(5) = tTp(5) =CTp(A)\{ci, C 2 ,cl}. 

By Proposition 14.II ii. the Schur complement S satisfies Assumptions (i) to (iv) 
of [EL041 Theorem 2.1] on each interval (cg-i, eg), I = 1, 2,..., L + 1. Hence the 
variational characterizations (I4.16P follow from [EL04I (2.9] applied on (cf_i,Cf), 
respectively; note that the values — c» and +c», respectively, of the functional p 
in [EL04I (2.3)] can be replaced by the end-points c^_i and c^, respectively of the 
interval (c^_i,Cf). 

The claims for the index shifts ng, in (cf_i,c^) follow from Proposition 14.II iii. 
Claim i) is the case £ = L -|- 1 in Proposition 14.11 iii). Claim ii) was proved in 
(I4.13P ; the proof of claim iii) is analogous to the proof of Theorem 13.91 iii). □ 

Remark 4.7. Define the monic operator polynomial 

(4.17) 7^(A) ;= (—1)^“^ ^ J^(c^ — A)'j5(A), domP(A) = dom5(A) = dom A, 

of degree L-1-1. One can show that V is weakly hyperbolic (comp. e.g. |La.u73inWf^ . 
[MarSSI § 31] for operator polynomials with bounded operator values); in particular, 
for every u G dom A, the function A i—^ (P(A)u, u) has L+1 real zeros which coincide 
with the numbers pi{u) < P 2 {u) < .. -PL+iiu) defined in (14.151) . 

The two-sided eigenvalue estimates for the one pole case in Subsection 13.31 were 
based on the solution formula for quadratic equations and do not readily generalize 
to the case of several poles. However, under stronger assumptions on the opera¬ 
tors Bg, we can estimate the gap (c^, cg + sg) in the spectrum of S in Proposition l4.1l 

Proposition 4.8. Let PL he a Hilbert space, H = 'H©'Hi©'H 2 ©- • ■(B'Hl with Hilbert 
spaces PL, PLi, PL 2 , ■■■, PLl- Let A he a self-adjoint operator in PL with compact 
resolvent and bounded from below, A > a for some a G M., Bg : PLg ^ PL bounded 
linear operators with the same closed range, ran.Bg =: PLq <G PL, i = 1,2,L, and 
Cl, C 2 , ■ ■ ■ ,cl G M with Cl < C 2 < • • • < Cl. Then there exist 7 ^ > 0 such that 

(4.18) jgWuof <\\B*gufg<j+\\uof, u = uo + u'oGno(Bn^, £= 1 , 2 ,...,L. 
Furthermore, for i G {1, 2,..., Lj the following hold: 

i) If a> c^+i and rjg is the unique zero of the function 

^ _ L + 

E AeR\{ci,C 2 ,...,CL}, 

^ CL - A CL - A 

in {cg,cg+i), then {cg,r]g) C p{S). 

ii) If Ho = H and 77 } is the unique zero of the function A 1 —>■ a — A — '0^(A) in 
{cg,cg+i), then {cg,r]g) C p{S). 

Proof. Let £ G {1,2,...,L}. By definition of Ho, the operator Bg : Hg ^ Ho is 
surjective. If we decompose H = Ho® Hq with Ho = ran Bg, Hq = {ran Bg)-^ = 
VerBf, then Bfu = Bfuo for u = uo + Ug G Ho ® Hq and Bf\no ■ Ho ^ Hg is 
injective. Further, since ranBg = Ho is closed by assumption, so is ran{Bf\'Ho ) by 
the Closed Range theorem. Hence, by the Closed Graph theorem, (RUho) ^ is 
bounded on ran{Bf\'ggg), and the existence of the numbers 7 ^ in (I4.18P follows. 

For u G domA, we consider the function 

Pu{X) ■■= {S{\)u, u) = ((A - X)u, u) - E 

. C£ — A 
n _1 ^ 


A G ]R\{ci,C2,...,Cl}, 
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where || • denotes the norm in T-Lt. The claims i) and ii) follow if we show that in 
a right neighbourhood of q all functions u G dom^l, are uniformly positive. 

For A G (cf,c^+i) we can estimate 






\Biu\\l 

k = l ktt +1 

L + 


>((,4-AK.)- y: 


7fe 


Cfe — A Cfc — A 

k=l k=e+i 


koll^ 


= {{A - X)u,u) - V'KA)||mo|P- 


Since is strictly increasing in (c£,q+i) from —oo to +oo, it has precisely one 
zero 77 ^ G (c£, c^+i) and it is negative on {ci,r]g). This implies that, in case i), 

(p„(A) > (a- c^+i)||m|P, Ag(c^, 77 ^), 

and, in case ii), with rj'^ as defined there, 

(a - A - V'^(A))||it|p > 0 Ag(c^, 77 ^). □ 


5. Application to photonic crystals 


In this section, we consider electromagnetic waves (E, H) propagating in a non¬ 
magnetic medium with relative permittivity e. The permittivity (or dielectric) 
function e depends on the spatial coordinates xi, X 2 as well as on the frequency w, 
but not on the spatial coordinate 0 : 3 . The electromagnetic wave {E,H) is then 
decomposed into transverse electric (TE) polarized waves {Ei, E 2 ,0,0,0, H^) and 
transverse magnetic (TM) polarized waves {0,0, E 3 , Hi, E[ 2 , 0 ) |DL901 Chapter 1], 
[FK96) . This decomposition reduces Maxwell’s equations to scalar equations for 
H 3 and i? 3 , respectively. 

The discussion below focuses on the TM case. We will apply the operator the¬ 
oretic results developed in the preceding sections to a physically relevant spectral 
problem with periodic permittivity function. This operator formulation of the pho¬ 
tonics crystal problem with A-dependent material properties is used as a base for 
numerical approximations in Section |6l 


5.1. The physical problem. Let A := 

operator in The equation 


dxl 


92 

-I- denote the formal Laplace 
0X2 


(5.1) CE 3 :=—AE 3 — u}'^e{x,uj)E 3 = 0, x := (xi, 0 : 2 ) G 


models transverse magnetic (TM) polarized waves with frequency uj. The permit¬ 
tivity function e will only depend on G P C C and we therefore define A := 

Let F denote the lattice and fl := (0,1]^ the unit cell of the lattice F. The dual 
lattice to F is F* := 27rZ^ and we define the fundamental domain (the Brillouin 
zone) of the dual lattice F* as the set /C := (—7r,7r]^. In this paper, we consider 
Bloch solutions of (EB, i.e. non-zero solutions of the form i? 3 (x) = e'^^’^'^u{x), 
X G where m is a F-periodic function, k € 1C, and u is a function of the variable 
X on the torus jl? [Kuc931 p. 104]. Here (•, •) and | • | denote the Euclidean 

inner product and norm in respectively. Further, let V denote the gradient 
with respect to the space variable x G R^, (f, V) the directional derivative in the 
direction v G R^, v ^ 0, and iJ®(T^) the Sobolev space of order 2 associated with 
L^(T^). Note that functions u G iL®(T^), s > 0, can be characterized in terms of 
their Fourier series with coefficients u{n) G C, n G Z^, 
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M S ^ (1 + |27rn|)^®|u(n)p < oo. 

For fixed k G K. = (—7r,7r]^ the shifted Laplace operator : L^(T^) —>■ L^(T^) 
is defined as 

(5.2) Afc := (V + ifc.V + i/c) = A + 2i(A:,V) - |fc|2, dom(Afc) = iL2(T2). 

The operator A^ has compact resolvent and a{—Ak) = {|27rn + /cp : n G 
[Kiic981 p. 161-164]. Since V(e'^^’^^u(j;)) = -|-ifc)u(a;), the Bloch solutions 

E 3 are formally determined by the solutions of Tk{)^)u = 0 with 

(5.3) rfc(A) :=-Afc-Ae(-,A), fc S/C = (-tt,T rj^, A G P C C, 

over the torus T^. The solutions u and eigenvalues A depend on k G K. but we 
will not write this out explicitly. A given A = is called a band gap frequency if 
equation (EH), independent of the parameter A: G /C, has no non-zero Bloch solution. 

General analytic properties of the permittivity function e(a;, •) are discussed in 
[Ces961 Chapter 1], while |EnglO| considers basic spectral properties of bounded 
holomorphic operator functions with applications to photonic crystals. 

Here, we apply the theory developed in the previous sections to a rational oper¬ 
ator function with periodic permittivity e. Let H = HiU • • • UHm, M gN, denote a 
partitioning of H = (0,1]^ and let denote the characteristic function (indicator 
function) of the subsets C (0,1]^, m = 1,2,..., M. In most applications the 
function e is piecewise constant in x, 

M 

(5.4) e(a;,A) := ^ em(A)xn,„(a;), xGfl,XG'DcC. 

m —1 

The A-independent case Cm ■= a-m leads to the extensively studied linear problem 
[FK961 rKuc93| . However, a common model for solid materials with dispersion is 
the Lorentz model 

Lrr^ , 

(5.5) e^(A) := + V ^ , A G C \ {c^,; : ^ = 1, 2,..., L„}, 

with Lm G N and positive constants Um, bm,i, and Cm,i |Ces96) . Here the permit¬ 
tivity function in (15.41) takes the form 

M M Lm , 

(5-6) e(-, A) := ^ amXOm (’) + ^ (•) 

with M G {1, 2,..., M} and Em '■= am, m = M + 1, M + 2,. ■ ■, M. By a polynomial 
long division of Ae(-,A) with the Lorentz model (15.6L we obtain 

M M Lm M Lm , 

(5.7) Ae(-, A) = A ^ amXn„(-) - ^ ^^m.£Xn„(-) + ^ ^ ^’^ XOmi') 

m —1 m —1 i —1 m —1 £—1 ’ 

for A € C \ {cm/ : ^ = 1, 2,..., m = 1, 2,..., M}. The operator 

M 

(5.8) W : L\T^) ^ L\T^), IT ^ a^xa„., 

m —1 

is bounded, self-adjoint, and bijective since >0. Moreover, 'Ylm=i Xoim 
is the identity operator in L^(T^) and {W*)~^ = W~^ = Ylm=i ^rn^XOm- 
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For A e C \ {Cm,£ ■■ £ = 1,2, ... ,Lm 
Sk{X) : L^(T^) ^ L^T^) by 

(^■9) 5fc(A) := W-^TkiX)W-^ = A 


,m = 1 , 2 ,..., M}, we define the operator 

M Lm 


X A) 


— 1 


-Xn„ 


m —1 ^—1 


where 

(5.10) Ak:=A]°^+A<^^\ 4°^:=-VF-^ AfelF"^, 

with domains 


M Lm 


dom5fc(A) = domA/c = dornA^^^ = domA^ = 

The operator is self-adjoint with discrete spectrum, while A^^^ is self-adjoint 
and bounded. Hence, A^ is self-adjoint and the spectrum a{Ak) is discrete as well. 


5.2. Block operator matrix formulation. In this subsection, we consider the 
minimal linearization of the rational operator function in (15.91) . Here the charac¬ 
teristic function xOm be viewed as a multiplication operator = Xn„- between 
L^(T^) and its range 

(5.11) ranPm = {u € L^(T^) : u\q\^^ = 0} . 

Then Sk can be written in the form 


M 

(5.12) Sk{X)=Ak-X-Y,J2^--xicm,e-X)-^B*^^„ 

m—1 1—1 



with Ak as in ()5.10l) . The theory in the previous sections assumes that the poles of 
Sk are disjoint and ordered increasingly. Therefore, we denote by L the number of 
disjoint poles of Sk , 

L .— ^ • m — 1,..., M, £ — 1,..., j- , 

and we denote these disjoint poles by ci < C 2 < • • • < cz,. For i = l,2,...,L, the set 
e {1,... ,M} : = t{m) e {1,2 ,.. .,Lm},Cm,t = Cij =: {mi, ■ ■ ■,^Ani) : 

consists of all indices m such that the permittivity in (15.51) has a pole at c^. Let 


P:=L2(t2), P :=Pi0P2---©Hl, 4 := 0 ranP^, 

meMi 


and define the operators L^(T^) and C : H ^ H hy 


(5.13) 



(5.14) 

C:= diag C 2 l^^ ... 


where 



(5.15) 



Then Sk is the first Schur complement of the block operator matrix Ak in P := 
L^(T 2 ) ©P given by 


Ak 


[Ak B\ 
\B* C J ^ 


domAlfe =domHfc © TL. 


(5.16) 
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Hence, the entries in the block operator matrix Ak are 



/ Ak 

)meMi 

{,^rn,l^{rn))rn^M2 




ciIfL^ 

0 

0 

Ak = 


0 


0 



0 

0 



Example. Let H = f2iUi72Uri3, i.e. M = 3, and consider a Lorentz model 
e(-,A) ■.=alxn^{■) + a2XnA-) + 


(5.17) 


^ 1.1 


ci.i — A 


Xni(‘) + 


-' 1,2 


Cl. 2 — A 


XOi(-) 


:^ 2.1 


C2,l — A 


XnA-) 


with Cl,2 < ci,i = C2,i. Here M = 2, L = 2, ci = ci,2, C2 = ci,i = C2,i, Mi = {1}, 
t{l) = 2, M 2 = {1,2}, £2(1) = £2(2) = = Hi,2, and 

= (Hi,iH 2 ,i). Hence, the operator function Sk is the first Schur 
complement of the block operator matrix 


Ak = 


/ Ak 

Hi,2 

iHi,i 

H2,1 

^1,2 

Cl 

1 0 

0 

Hr,i 

0 

1 C2 

0 

V ^2.1 

0 

1 0 

C2 


\ 


domAlfc = domAfc © (Hi ©7^2), 


in H© (Hi ©H 2) = L2 (T2) © (ran xoi © (ran xoi ©ran xoa))• 


5.3. Spectral properties. In this subsection, we analyze the accumulation of 
eigenvalues at the poles of the permittivity function by means of the results of 
Section |3l 

Let {i/j(Al^°^)}°^^, denote the eigenvalues of AM in (I5.10L 


respectively. Since am,e., bm,e > 0 and J2m=i we have 

M Lm 

0 < = sup 'Y' (— bjn,ePmU, u) 

||u||=i - ^ > 


Mil. \ —^ 

< max — > bm, 

m=l \ a. 


e=i 

M 


e=i 

Lm 


I ) sup 


'Y^{PmU,u) 


M / 1 ^—A 

< max 1 — > bm, 


m=l \ a„ 


e 


e=i 


Hence, by means of the classical min-max variational principle [R.S781 Theorem 
Xlll.l], the eigenvalues of the operator Ak defined in (15.101) can be estimated by 


(5.18) 


\k\^ < Vj{AY) < VjiAk) < VjiAY) + nmx — bm,i ■ 

Tfl —i 1 (Xrn . I 




In order to bound the norm of B defined in (15.131) . we use the inequality 

' ^m,£{m) bm,^(rn) 


ll(H,„.<!qm))mGMj| < Hiax 

' ' meMi 


E 11^- 


m^Mi 
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and \\Prn\\ < 1 to obtain 


(5.19) 


|B|| < V ||(S„,£.(m))meMj| < V niax 

m^Mi 


'^m.Um 1 ^m.Um) 


Clearly, dim 'Hm = oo and the operators : ranP^ —t are injective. Take 

i € {1, 2,... ,L} arbitrary and assume Ui = (itm)meMi € kerP^, i.e. 


m^Mi m^Mi 

Then P^Um = 0 since the coefficients of P^Um are positive, the support of P^Um 
is contained in VLrm and the subsets are pairwise disjoint. Hence Um = 0 for all 
m S Mi, i.e. Ui = 0, which implies that Bi, z = 1, 2,..., L, is injective. 

Now it follows from Propositions 14.31 and 14.41 that all poles of the permittivity 
function e in CT . i.e. all points Ci < C 2 < ■ ■ ■ < cl, are accumulation points of 
eigenvalues of Ak- Moreover, Proposition 14.11 implies that, for fixed k, there are 
gaps in the spectrum of Ak to the right of all a. Note that only depends on 
k and the dependence is analytic. 

In general, we will not have a band gap above the poles for all material models 
(15.61) since the width of the gaps may tend to 0 if fc varies, i.e. we may have 
A^+i,i(fc) ce, k ^ ko, for some fco S (—tt, tt)^ and i £ {1,2,..., L}. An example 
for this is given in Section |6] where Figure [5] shows that numerically we do not find 
a band gap above Ci = Ci_i. 

However, Proposition 14.81 provides a concrete estimate for a gap above a pole ci 
in the case when Ak > c^_|_i and all operators B^ i, £ = 1,2,..., L, have the same 
closed range. Since the range of Bm,e equals the range of P^, which is 

(5.20) Ta,nBjn,i = ranP,{j = [u £ L^(T^) : = O}"'", 

operators Bm,i, Bn/ with m n do not have the same range. Thus Proposition l4.8l 
only applies to permittivity functions of the form 

M P h 

(5-21) e(-. A) := ^ a^xn„ (•) + ^(•) 

m=l t.=\ 

and to poles c\/ such that Ak > Ci/+i. Since e{-, A) and hence the operators Bm/, 
do not depend on k. Proposition |T8]i) shows that if Ak > k £ JC = (—7r,7r]^, 

then there exists an r/ > ci, independent of k, such that 

(5.22) (q, rje) C p{Ak), k £ 1C = (-tt, irf, 

which means that in this case our abstract results guarantee a band gap above c^. 


-- —PmUm- 


Explicit bounds on the band gap above eg,. In this subsection, we consider a case 
where an estimate of the band gap in (j5.22l) can be derived explicitly. Assume the 
permittivity function has the form (j5.21|) with M = 2, Li = 2, and set ci := ciy < 
Cl,2 =: C2. 

Let u = uo u'q £'Hi ® 'Hi - By the definition of 
(5.23) = 

Oi 
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which implies that the constants in Proposition 14.81 are 
(5.24) 


_ biiCii 

7i = 


Oi 


+ &1.2Cl,2 

72 = ' 


Oi 


According to Proposition 14.81 we let 


(5.25) 


•02 (A) = 


7i 


7.+ 


ci.i — A ci ^2 — A ’ 

It is easy to see that the unique zero iji of ip 2 in the interval (ciq,ci^ 2 ) is given by 

Cl,1^1,1 


(5.26) 


m = ci,i + 


A G R \ {ciq, ci_2}- 
the interval (i 
■(ci,2 — Cl,l). 


Cl,2&l,2 + Cl,l6l,l 

Hence, for ci ,2 such that Ak > ci ,2 for all fc G /C = (—7r,7r]^, Proposition 14.81 iil 
implies that we have a band gap above cip, namely 

Cl,1^1,1 


(5.27) cpi,cpi + 


(ci, 2 -cip) Cp(7lfc), fc G /C = (-7r,7r]2. 


Cl,2&l,2 + Cl,l6pi 

The accumulation of eigenvalues at the poles of the permittivity function from the 
left and two-sided estimates for them will be discussed and illustrated in Section [51 

5.4. Space independent permittivity. The case e{x, A) = e(A), where the eigen¬ 
values of Ak can be calculated explicitly, gives additional insight into what we can 
expect from the numerical calculations in Section 16.31 In particular, it shows that 
all the abstract two-sided eigenvalue estimates in (13.251) are optimal. 

Let a, b, and c be positive constants, 

(5.28) e(A):=a-b—A gC\{c}. 

c — A 

Here M = M = 1, Pi = Ji 2 (T 2 ), and W = a/i 2 (T 2 ) since H = Hi. Hence Ak = 
— iAfc -I- ^ and B = Bi = y^Ji 2 (T 2 ) so that the estimate in (I5.19|) is, in fact, 
an equality, ||i?|| = and we know the eigenvalues of Ak and of Ak explicitly. 

With the notation {vj{—Ak)}'^i = {127rn-|-fcP : n G I?} introduced earlier for the 
eigenvalues of — in increasing order, we have 


(5.29) Wj{Ak)}%i = { 


Vj{-Ak) + b\ 




and the two sequences of eigenvalues of Ak accumulating at c and at oo are given by 
{Ai.i}r=i = {i(^i(Alfc)+c)-^i(:.i(A,)-c)V^|“^, 

{A2,i}^i = ^-^(tyj{Ak)+c^ + )J -(^i^j{Ak)-cj +~| ■ 


j=^ 

I'l 


The order of convergence of Aij- —?> c and of X 2 J —>■ 00 is 0{i^j{Ak)~^) as expected 
from Remark 13.171 

To compare with the two-sided eigenvalue estimates in Theorem 13.131 we need 
to determine the index shift K 2 from Proposition l3. 61 which is defined as the number 
of negative eigenvalues of 5^ (c -|- (5) for any <5 > 0 such that c < i5 -I- c < A 2,1 . It is easy 
to see that the eigenvalues fj,j of Sk{X) = Ak — X — for A = c -f A are given by 

r c6 1 1 °° 
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which implies that all eigenvalues are positive if <5 > 0 is chosen small enough and 
hence K 2 = 0. 

Since B — Bi = y^Jj;, 2 (T 2 ), we have mma{BB*) = ||i?|p = and the two- 
sided estimates (13.251) all become equalities, i.e. 

(5.31) = Ap, = A^^., X^^ = A2., = A^^^., 

which proves that they are optimal. 


6. Galerkin approximations 


In this section, we consider Galerkin finite element approximations of eigenpairs 
(A, It) of the operator function (15.31) with permittivity function (15.61) or, equivalently, 
eigenpairs of the Schur complement (15.91) . Since, in Section lOl 16.21 below, we use 
results that were derived for bounded forms, the forms associated with the operators 
(15.31) will be considered in i?^(T^) rather than in L^(T^). Moreover, here it is more 
convenient to use the operator W in (15.81) in a new inner product and work with 
7/c instead of Sk ■ 

In the space iJ^(T^) x i?^(T^), for fixed k G 1C = (—7r,7r]^, we consider the 
bounded sesquilinear forms 


4°^ [w, u] := / (V -b ik)u ■ (V -b ik)v dx, 

Jn 

(6.1) ik[u,v] := -b 



M L„i 

:= ^ y^bm,ibm[u,v], 

m—1 ^—1 


M 

{u,v)w ■= ^ am,bm[u,v], 

m—1 


and denote by ||u||u; = \/(u, u)u, the corresponding weighted L 2 -norm. Then the 
eigenvalues of 7/c are determined by the following variational problem: 


Find u G i7^(T^)\{0} and A G R\{cm,£ : to = 1, 2,..., M, I = 1,2,, L^} such 
that for all v G 71^ (T^) 


M 

(6.2) t/c(A)[M, u] := ik[u,v] - X{u,v)y, - V] ^ Cm./; = q. 

From the preceding sections, we know that the eigenvalues of Sk, and hence of Tk, 
are isolated and of finite multiplicity and hence so will be the values A for which 
the above problem admits a solution (comp. e.g. |BS081 p. 131]. 

Let denote an arbitrary TV-dimensional subspace of 77^ (T^). Then the dis¬ 
crete (Galerkin) formulation of problem (16.21) is to seek the N eigenpairs (Xj,Uj) G 
(M\{cm,f : TO = 1, 2,..., M, 1 = 1,2,..., L^}) x 'Hjy such that for all f G 77^ 

(6.3) t/c(Aj)[Mj,z;] = 0, j = 1,2,...,7V. 

For a sequence of subspaces C 77^ (T^), the corresponding sequences of eigen¬ 
pairs {Xj,Uj), j = 1,2,..., N, are viewed as approximations to the true eigenpairs 
{Xj,Uj) of the spectral problem Sk{X)u = 0 (comp. e.g. [BOSQ] !. 

According to Theorem 14.61 all eigenvalues of Sk, and hence of Tk, can be deter¬ 
mined by means of the variational principles (I3.20L (13.211) . For a material model 
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dSll) with M = 2, M = 1, and Li = 1, the functionals pi ,2 in (13.1911 are given by 


(6.4) pi,2(w) = 


1 / tk[u,u] 


+ ci ,1 T 


1 / tk[u,u] 


- Cl,l 


bi^ici^ibi[u,u] 


for u G Note that, in analogy to Remark 13.101 the variational principles 

(I3.20|) . (j3.21ll hold with domAfc = replaced by domtfc = 

In the same way as for the eigenvalues of Sk , we divide the Galerkin eigenvalues 
into two groups. By 

(6.5) Aip < Ai^2 < • • • < ^i.Ni < cin 


we denote the Ni (< N) eigenvalues of (16.31) below cip and by 

(6.6) Cip < A2,1 < A2,2 < • • • < ^2,N2 

the N 2 (= N — Ni) eigenvalues above cip; the corresponding eigenvectors are 

denoted by (m 2 j)^i Then the variational principle (I3.20p for 

together with Remark 13.101 implies that 

(6.7) Ai j = min max pi(m) < max ^ Pi(u) < Ai j = 1, 2,..., Ai, 

■uespan{ui,...,uj-} 

dimZ:=j 

and, similarly, the variational principle (13.211) implies \ 2 ,j < A 2 J , j = 1, 2,..., fV 2 . 
So, in analogy to the classical min-max (or Rayleigh-Ritz) variational principle, we 
have the chain of inequalities 

(6.8) Ai j- < Ai j < cip < A2,fe < A2,fc, j = 1,2,..., Ni, k = 1,2,... ,N — Ni, 

for the eigenvalues of 7fc and the Galerkin eigenvalues to the left and to the right 
of the pole ci 1 . 

Regarding the eigenfunctions, let G K x i7^(T^), j G N, be a sequence 

of eigenpairs of Sk, here given by 

(6.9) 5fe(A)=Alfc-A-Ri,i(ci.i-A)-iRri, Bl,:= 

’ ’ V ai 

such that Aij ^ cip, j —>■ 00 , and ||uij|l = 1, j G N. Then the corresponding 
eigenvectors are of the form (uij uijY with uij = —(cip — X)~^B*uij, j G N, 
and {(ui_j is a singular sequence of Ak at cip, i.e. 

||(-4fe —ci,i)(uij-uij-)*|P= ||(Afc —cip)uij- + iluij| 1 ^ 2 C]p 2 ) + ||7?*uij||^—>-0, j —>■ 00 . 
This implies that 

(6.10) ||i?*uij|| ~ = 6c j |uij(x)|^dx —>■ 0, J —>■ 00 . 

iOi 

Thus we may expect that, for large j, the L 2 -norm of the eigenfunctions uij of 
Sk on fli is very small. We mention that while the order of convergence for the 
eigenvalues Aij- ^ cip is 0{i'j{Ak)~^) according to Remark 13.171 the order of 
convergence \\B*uij\\^ —^ 0 will depend on the geometry and other parameters. 
Furthermore, since (ci,i,cip +e) C p{Ak) for some e > 0, we may not expect that 
the I |il*U 2 j I |z, 2 (f 2 i) is small for eigenvectors U 2 ,j of Sk corresponding to eigenvalues 
X 2 J above ci.i. 

Figure [3] illustrates this by numerical computations for a particular example 
where cip = 12.7367. The constants in the material model (EH) are M = 2, 
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M = 1, Li = 1, and is a disk of radius r = 0.475. Numerically the convergence 
—>■ 0 is faster than (ci,i — —>■ oo in this case, which can be seen 

in the third plot from the left. The numerical method used here will be outlined in 
the next two sections. 



Figure 3. fli disk of radius r = 0.475 (black circle) and pole ci,i = 

12.7367. The absolute value of the approximate eigenvectors u cor¬ 
responding to (counting from the left) the eigenvalues Ai,i = 1.4012, 

Ai ,2 = 11.6218, Ai ,3 = 12.613, and A 2 ,i = 14.3949. Blue corresponds to 
|u(a:)| close to zero and red to large function values. 

In the following numerical computations with space-dependent permittivity func¬ 
tion e we will approximate the analytic bounds on the eigenvalues of Ak established 
in Theorem 13.131 The estimates (I3.28|) , p.29p only involve the eigenvalues of the 
positive densely defined unbounded self-adjoint operator Ak, which has compact 
resolvent; note that we already know Tnh\a{BB*) = 0 and ||S||, or can estimate 
the latter. The spectrum of the operator Ak can be obtained from the following 
variationally posed eigenvalue problem: 

Find u G i7^(T^)\{0} and u G R such that for all v G H^(T^) 

(6.11) tfe[u, w] = u(u,u)w. 

The Galerkin approximations i/j of Vj{Ak) G a{Ak), j = 1,2,..., N, with ui < zAj < 

■ • • < and the corresponding eigenvectors are eigenpairs (u, u) G R x 'H^\{0} 
such that for all u G 

(6.12) tfe[w,u] = u(rt,F)w. 

The classical min-max variational principle |RS781 Theorem XIII. 1] applies and 
shows that the Galerkin eigenvalues are upper bounds for true eigenvalues |RS781 
Theorem XIII.3], 

(6.13) u,(Alfe)<u„ j = l,2,...,N. 

Therefore, by Remark [3311 since the bounds A^^ and are increasing func¬ 

tions of Vj{Ak), the eigenvalues satisfy the estimates 

(6.14) Xi,,<Xl,<~Xl,, 

(6-15) A 2 J-I-K 2 ^ + 

where Ai ^ and X 2 j are defined by (I3.28P and (13.291) . respectively, and where A^^-, 
are obtained from X^ j by replacing Vj{A) by its upper bound Vj. 

6.1. The conforming finite element method (FEM). The conforming Galerkin 
FEM is characterized by a family of triangulations Qh of ft with quadrilateral el¬ 
ements Q. Let Fq be a bijective mapping of a chosen reference square Qref onto 
an element Q and let denote the space of polynomials on R^ of degree < p, i.e. 
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PP := span 0 < n, m < p}. Then the space of polynomials on of degree 

< p that are piecewise polynomials on the triangulation of is defined as 

(6.16) VI := {u e : u\q o Fq S P^}; 

here u\q denotes the restriction of m to Q- The topology of the torus is imposed by 
mapping the parallelogram edges situated on one boundary on the corresponding 
parallelogram edge on the opposite side. From the Cauchy-Schwarz inequality it 
follows that for t S [ 0 , 1 ] and n S 1 ? the inequality |27rn+fcp > (27r)^|np(l — t) + |fcp 
holds when jfcj < tir. Hence, for jfcj < tir with t G [0,1], the inequality 

(6.17) ife[M,u] > / |(V + iA:)Mpdx> (l-t)||Vu|p + |fcn|u||2, uGH\T^), 

Jn 

follows. Thus, for k 0, is a norm equivalent to the standard norm on 

and ife)!!,!;] can be used as an inner product. As above, let 
denote the eigenvalues of Ak ordered increasingly and let be a system of 

eigenvectors that are orthonormal with respect to the inner product tfc[-,-], i.e. 
Ifc [Hz; 1 bj ^ 

The corresponding finite dimensional problem is (|6.12l) with the special choice 
= Vj^, i.e. to seek u G 7^^\{0} and u G M such that for all v G V^ (16.121) holds. 
Let 0 < ui < U 2 < • • • < denote the eigenvalues of (16.121) with corresponding 
eigenvectors ui,U 2 , ■ ■ ■, um orthornormal with respect to ife[', •], i.e. ik^i, Uj] = 6ij, 
i,j G N. Moreover, let P denote the orthogonal projection of i7^(T^) onto V^- 
Then, the well-known estimate [Wei741IKO06] 

(6.18) 

i=l 

holds which, under the assumption that X]i=i IK* “ ^ tie written in 

the equivalent form 

Ifiiql Uj - UjjAk) ELi \Wr - Pu^\\l 

Note that ()6.19l) depends on the approximation errors for all eigenfunctions ui,U 2 , 

... ,Uj. It is possible to derive error estimates that depend mainly on the approx¬ 
imation error for Uj [KO06] but this will not be important in our setting since all 
eigenvectors are (piecewise) analytic. In the following we will express the conver¬ 
gence rates in terms of the number of degrees of freedom N. If Uj is analytic, the 
convergence rate for the p-version of the finite element method is exponential, 

( 6 . 20 ) 

for some positive constants C and 7 |BG 88 ) . Moreover, the estimate (16.201) holds 
when the eigenfunctions are analytic in each subdomain up to the interfaces and 
the interfaces are exactly resolved by curvilinear cells |BG 88 | . 

In the following, we discretize the operator Ak and the block operator matrix 
(14.21) using the implementation available in the software package Goncepts [FL02j . 
The shape functions are based on Jacobi polynomials and a blending technique 
is applied to construct element mappings. Goncepts uses curvilinear quadrilateral 
cells that resolve the curved analytic material interfaces in the numerical examples. 
Further implementation details can be found in |FL02l ISK09j . 

Let { 01 , (/ 12 ,..., be a basis of V^. For u G V^, the approximate eigenfunc¬ 
tions u are of the form 
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N 

( 6 . 21 ) u = 

_ _ 1=1 _ 

The elements in the matrices Tk,W, and are 

( 6 . 22 ) {Tk)ij = ik[(l>j,4>i], {W)ij = {Bm)ij = bm['Pj,4>i], 

where the sesquilinear forms are defined by dSI]). The finite element approximation 
of the rational eigenvalue problem (16.21) then is 

(6.23) 7fc(A)a; = 0, x = (ai, a 2 , • ■ •, a^)*, 
where, for A S C \ {cm,i : m = 1,2,..., M, £=1,2,..., Lrn}, 

_ ^ _ M Lm , 

(6.24) Tk{X) ■.= Tk-\W-J2J2 

The eigenvalues of Ak := are used to approximate the eigenvalues 

of Ak and hence the bounds on the eigenvalues of Ak in (13.281) . (13.291) . The rational 
eigenvalue problem Sk{X)x := W~^^'^Tk{X)W~^^‘^x = 0 with M = 1 can be written 
in the form (14.21) and we will calculate the approximate eigenvalues using this block 
matrix form. 

Remark 6.1. Note that Bm is not block diagonal for conforming methods and 
the method outlined in this section can only be used when M = 1. However, 
the approximation spaces for discontinuous Galerkin (DG) methods are localised 
in each element and Bm is then block diagonal. A DG method that can handle 
M > 1 is outlined in Section 

The variational characterization of the eigenvalues (IgTl) shows that (ignoring 
rounding errors) the approximate eigenvalues are upper bounds on the eigenvalues 
of Ak (see (16.81) '). 

Given two closed subspaces Hi and H 2 of a Hilbert space H the proximity of 
the spaces is measured in terms of the containment gap |Kat95[ Section IV.2.1] 

?(Hi,H2):=max{5(Hi,H2),(5(H2,Hi)}, sup inf 

For the gap 5{E^{X),E{X)) between the discrete and continuous eigenspaces E^{X) 
and E{X) of Ak, respectively, it is known that 6{Ef,^{X), E{X)) —?> 0 when /i—>-0 or 
p —> 00 [Osb75) . This implies several decisive properties including non-pollution of 
the spectrum [DP04) . Hence, without risk of spectral pollution, we can calculate 
accurate approximations of the estimates Af ^ in (13.261) and X^j in (]3.27p . We men¬ 
tion that a complete convergence theory for Ak will be studied in a forthcoming 
paper; a block operator formulation of a related problem was considered in |Engl4| . 
Since all eigenfunctions are analytic in each subdomain up to the interfaces we ex¬ 
pect exponential convergence. The presented continuous finite element method will 
be used in Section [Q to compute the eigenvalues but the discontinuous Galerkin 
method described in Section [6.21 below is a strong alternative. 

6.2. The symmetric interior penalty method. In the following, we discretize 
(lO) with a discontinuous finite element method called the symmetric interior 
penalty method (SIP) [ABP06| . More specifically, we use the p-version of SIP 
in this paper since in the numerical examples all eigenvectors are analytic in each 
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subdomain and the subdomains have analytic interfaces. Discontinuous Galerkin 
methods, such as SIP, are interesting partly because they are more flexible in the 
choice of basis functions and in the mesh design |Gial2) . Moreover, the mass matrix 
is block diagonal and can therefore be inverted at low computational cost. This 
is frequently used to obtain explicit time integration, but a block diagonal mass 
matrix is also an advantage for the solution of our non-linear eigenvalue problem. 

Let Th denote the triangulation of VL, let Ft be a bijective mapping of a chosen 
reference triangle TJ-ei onto an element T, and let denote the space of polynomials 
on of degree < p that are piecewise constant on the triangulation of D, 


(6.25) 


Vl := {u e : u\t o Ft £ P^}; 


here u\j. denotes the restriction of u to T. Gonsider two adjacent triangles r+, T_ S 
Th with outward pointing normals ri_|_, fi- on the shared edge i9T+ fl i9T_. The sym¬ 
bols Vh and Vh- denote the elementwise (broken) gradient operator and divergence 
operator, respectively. The solutions of (j6.2l) are in iJ^(n). Hence, for all T £ Th 
the traces on the edges are well-defined. The averages {•} and jumps [•] of w and 
V/iW on T+ U T_ are then defined as 


(6.26) {w} := |(ui+-I-ui_), 

(6.27) {yhw}:=^{VhW+ + VhW_), 


[w] := ■u;-|-n+ -I- W-fi-, 


where w± and Th'UJ± denote the traces on dT±. Let £ denote the set of all edges 
e of 7/j. We will use the convention 


(6.28) 


Uhdx := X! / / Uhds ■='^ Uhds. 


TGTh' 




For Uh,Vh £ Vf we consider the sesquilinear forms 


:= J" (^(Vh + ik)uh, (Vh + ik)vh'^ dx + J" (/][uh], [vh]) ds 

-J ({(V/j-b iA:)M/i}, [u/j])-f (|{(V;i-t-iA:)?;/i}, ds, 

ik[uh,Vh] ■■=^°\uh,Vh] +t^'^'>[uh,Vh]. 

Let {(^ 1 ,(^ 2 ,.. ■, be a basis of V^. The finite element matrices and the corre¬ 
sponding rational matrix function Tk are then formed as in (I6.22l) - ()6.24l) . with 


(6.30) 


(Tk)ij — ffc 5 (Bm)ij — ■ 


The periodic boundary conditions are imposed weakly by identifying opposite sides 
of the unit cell and enforcing periodicity of the solution via the corresponding 
penalty terms in (16.2911 . To generate the mesh, we use the software package 
Emc^ [SH95j and for the discretization we build on the Matlab version [HW08j 
of NUDG++ (www.nudg.org). In order to preserve the high accuracy of the p-version 
of SIP it is essential to use curvilinear elements. A Gordon-Hall blending [GH73] 
is used to preserve the approximation properties of the basis. One advantage of 
SIP is the block diagonal mass matrix which directly splits into mass matrices for 
rn = 1,2,, M. Hence, Bm, m = 1,2,..., M, are block diagonal matrices. 
Let, for example, M = 2 and write Bi, B 2 in the form Bi = diag(ZIi,0) and 
B 2 = diag ( 0 , 112 ), respectively. The matrices Di £ and £>2 £ ]R" 2 xn 2 


are 


















34 


CHRISTIAN ENGSTROM, HEINZ HANGER, AND CHRISTIANE TRETTER 


positive definite, where ni, n 2 correspond to the number of basis functions sup¬ 
ported in domains r22. Let Dj =: L*Lj denote the Cholesky decomposition of 
Dj, j = 1,2. Then the Cholesky decomposition of the block diagonal matrix W is 


(6.31) 


W--=%L^ 


■ [ 0 


0 


\ 0 


0 

1/02^2 


Note that, for a given polynomial degree p, Dj, j = 1,2, are block diagonal with 
block sizes (p -I- l)(p -I- 2)/2 [HWOSl p. 171] and Lj therefore has only triangular 
blocks on the diagonal. Hence, the inverses of Lj and of can be computed at 
low costs. Let Inj, j = 1,2, denote the identity matrix on and define 

Sk{X) :=(%)-i^(A)(L^)-i =3; - 1 -H 2 .i(c 2 ,i - 1 


where A, := B^,^ := 0), H 2.1 

The function Sk is the first Schur complement of 



(6.32) 



/ Ak 

Bi,i B2 ,i \ 

— 

Bli 

ci,i 0 


\Bli 

0 C2,l / 


Hence, the eigenvalues can be computed from the 2(ni -|- 712 ) x 2(ni -|- 77 , 2 ) ma¬ 
trix Ak- For the operator Ak convergence in the gap of the eigenspaces is known 
[ABP061IGH12) . which implies non-pollution of the spectrum. However, a conver¬ 
gence analysis of the discretization of Ak is beyond the scope of the current paper. 
The main aim of our calculations in Section [631 is to illustrate the general theory in 
the previous sections and to show how our abstract results apply in concrete cases. 
Note that, for both finite element methods, convergence theory is known for the op¬ 
erator Ak, but the bounds (16.81) only apply to the conforming finite element method. 


6.3. Numerical examples. In the last twenty years physicists and engineers have 
studied d ispersive (A-depen dent) materials in periodic structures extensively 
[SSCH941 lHB,l+03l lHL,T+04l ISKE06] . In particular, polaritonic materials have 
received much interest, mainly because they exhibit a strong resonance at infrared 
frequencies [GMKLl^ . A common model for polaritonic materials is 

(6.33) e(A) :=e^ ("1 + ^ 4 ^ 

Contrary to ordinary non-dispersive structures, the physical lattice constant d (the 
physical unit cell is dx d) plays an important role. Let d = 21/im denote the lattice 
constant and the speed of light. The material constants of gallium arsenide 
(GaAs) in units of ^-kv ja then are [GMKLl^ . 

(6.34) £00 = 10.9, WL = 0.612, wt = 0.568. 

In all examples the standard IRA algorithm as provided by ARPACK [LSY98) was 
used to compute the eigenvalues. Moreover, in the one pole case residual inverse 
iteration [Neu85] based on the Rayleigh functionals (I3.19|) was used to verify the 
computations of Aip and A 2 ,i defined by (16.5p . (|6.6I) . 

6.4. One pole case. The conforming finite element method in Section lOl was used 
to study the polaritonic material model (16.331) when the constants in the material 
model (15.6p are M = 2, M = 1, and Li = 1. The material in Hi is polaritonic 
with constants (16.341) and H 2 is filled with air, e = 1. Hence, the constants in the 
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material model 

e(x, A) 
are 



^ 1,1 \ 
A ~ C14 J 


XnAx) +a2Xn^{x), 


a:Gr2=niun2, 


(6.35) ai = £ 00 , 02 = 1, 61,1 = eoo(a;L - Wt), ci,i=a;T 

with £ 00 , wl, and wt as in (16.341) . 

Table [T] shows that the approximations ui, V 2 to the first two eigenvalues i'i{Ak), 
V 2 {Ak) of Ak decrease for higher polynomial degrees p, which is expected from the 
classical min-max principle ()6.13|1 . The numerical calculations indicate that the 
convergence is exponential and that the estimate (16.1911 with (16.201) holds. Moreover, 
Table [T] suggests the same convergence behaviour for the approximations to the 
lowest eigenvalue Ai^ of the operator matrix Ak and to the first eigenvalue A 2.1 
above Ci^. This behaviour is expected from the variational principle (I3.20L (16.71) . 


Table 1. Polaritonic material model (16.331) . Eigenvalues of Ak and 
Ak when k = (vr,0), a\ — 10.9, a 2 = 1, 61,1 =22.3419, ci,i = 12.7367. The 
bold numbers show the digits in common with the solutions for p= 10 . 


p 

N 

Vl 

V2 

Ai,i 

A 2,1 

4 

320 

3.544840275 

5.168021445 

1.402494821 

14.39621821 

6 

720 

3.544596578 

5.166804859 

1.402239773 

14.39603122 

8 

1280 

3.544595950 

5.166797931 

1.402238958 

14.39602980 

10 

2000 

3.544595948 

5.166797908 

1.402238956 

14.39602979 


Table [2] shows the numerical approximations of the bounds (I3.27I) - (I3.29I) for a 
few eigenvalues when dim'H^ = 7320 and dimUfj = 3504. Recall that N{Ak,c) 
denotes the number of eigenvalues of Ak in less than or equal to c. In the 
example we have N{Ak,cip) = 4. Hence, for the matrix problem the condition 
dim'Hjy — diiR'Hfj > N{Ak,ci^i) in Remark 13.121 ensures that the shift K 2 for 
the finite dimensional problem is computed correctly. In this example, numerical 
calculations show that 5(c)''does not have any negative eigenvalues, which implies 
that the index shift K 2 is zero. 

Note that the chosen material parameters (16.341) give ||H|| = y^ 6 i,iCi 4 /ai ~ 5.1, 
which is relatively large; in this case we cannot expect the bounds (13.251) to be tight. 

Table 2. Polaritonic material model (16.331) . Bounds on Aij', A 2 .j for 


Sir disk of radius r 

= 0.475, 

k = 

(t,0). 

N{Ak, ci,i 

) = 4, 

0,2 = 1, Cl,l = 

12.73668500, and here 61,1 = 

= 22.3419, 

K2 = 0 . 

n 

^ 1,7 

Ai.j 


'^ 2,7 

A2.i 


1 

1.25 

1.39 

3.53 

Cl.l 

14.38 

15.01 

2 

2.57 

3.09 

5.14 

Cl,l 

14.82 

15.31 

3 

5.93 

6.04 

9.76 

Cl,l 

16.44 

16.57 

4 

7.00 

7.30 

11.55 

Cl,l 

16.96 

17.29 

5 

9.42 

9.56 

Cl,l 

17.31 20.48 

20.63 


Table [3] shows the numerical approximations of the bounds (13.251) for a few 
eigenvalues when 61 j = 1 , which gives a less strong rational term || i3|| Ri 1 compared 
to the data used to compute the values in Tabled Notice that the bounds in the 
latter case are much tighter. The shift K 2 = 2 is computed numerically as above. 
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Table 3. Polaritonic material model (16.331) . Bounds on Aij', X 2 ,j for 
fli disk of radius r = 0.475, k = ( 7 r, 0 ), N{Ak,ci^i) = 4, ai = 10.9, 
0.2 = 1, ci,i = 12.73668500, and here foi,i = 1, fv 2 = 2. 


n 


^ 1,1 


"' 2 ,. 7+2 


"' 2 ,. 7+2 

1 

1.60 

1.61 

1.71 

Cl.l 

12.96 

12.97 

2 

3.54 

3.57 

3.67 

Cl.l 

13.05 

13.09 

3 

7.67 

7.68 

7.90 

15.47 

15.82 

15.85 

4 

9.47 

9.51 

9.83 

16.42 

16.67 

16.72 

5 

12.36 

12.38 

Cl,l 

20.07 

20.21 

20.23 


6.5. Multi-pole case. For simplicity we consider only the case of two poles, where 
the constants in the material model (j5.6|) are M = 2, M = 1, Li = 2, and fli is a 
disk of radius r = 0.2. The symmetric interior penalty method is used to discretize 
the block operator matrix (14.21) . A few eigenvalues Aij(fc), \ 2 ,j{k) are numerically 
calculated using a selection of vectors k along the line segments between the points 

(6.36) r = (0, 0), A = (tt, 0), M = (tt, tt). 

The triangular path formed by these points is called the boundary of the irreducible 
Brillouin zone [,T,TWMO^ . 



Figure 4. Lorentz models (|6.37|l with 2 poles: Band structure with a 
disk inclusion fli with r = 0.2. The dashed line marks the upper limit of 
(15.2711 . In both cases ai = 4, a 2 = 1, &i,i = 8 , ci,i = 1; above 61,2 = 10, 
Cl, 2 = 2, below foi ,2 = 2.53, ci ,2 = 5. 






















NON-LINEAR EIGENVALUE PROBLEMS AND APPLICATIONS TO PHOTONIC CRYSTALS 37 


Figure [4] shows a few eigenvalues for the model ai = 4, 02 = 1, &i,i = 10, 
ci,i = 1, 61,2 = 4, Cl ,2 = 2. The dashed line corresponds to the upper bound 
in (15.271) . Note that these bounds only require that Ak > ci ,2 holds. Hence, in 
this case we can guarantee a band gap using for instance the verified eigenvalue 
enclosures in [HPW09| to show that Ak > ci, 2 . This problem will, in general, be 
much less demanding than proving the band gap directly. 

Lastly, we consider a case of a permittivity function (15.611 where both materials 
are A-dependent 

(6.37) €{x,X):=(ai+ ) XOi(x) + (^ 02 + ) XOsa;Gf2 = HiUH 2 . 

V Cl,1-Ay V C2,i-Ay 

Here Proposition l4.8l does not apply since Hi,i and B 2.1 do not have the same range. 

Indeed, the claim of Proposition 14.81 does not seem to hold. Figure [S] shows that 
numerically we have a fc-dependent eigenvalue above ci,i that tends to ci,i if fc— 
for some values on fco G So, in this case, we cannot expect a band gap. 



Figure 5. Lorentz model (16.371) with 2 poles: Band structure with 
a disk inclusion fli with r = 0.2, ai = 1, a 2 = 4, 61,1 = 8 , ci,i = 1, 

62,1 = 12, C 2 ,i = 2. No band gap above ci,i = 1. 
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